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CLAUSIUS-MOSOTTI EQUATION FOR CERTAIN TYPES 
ANISOTROPIC 


ABSTRACT 


The Clausius-Mosotti equation derived for crystals whose particles are 
anisotropic and arranged any lattice provided that all the particles are the 
same type and provided that the principal axes their tensors polarizability 
are parallel each other. The formulas are special interest for artificial 
dielectrics used lens material for microwave frequencies. They explain 
experimental results. 


The Clausius-Mosotti equation (CME) establishes relation between the 
dielectric constant substance, the polarizability and the density 
its particles. Derivations for crystals cubic lattice and for amorphous 
substances are well known. 

van -Santen and Opechowski (6) have derived generalized equation for 
crystals containing particles different types provided their polarizabilities 
are independent the direction the field. Jonker and van Santen (2) 
showed that this equation was not sufficient explain the essential properties 
barium titanate. Slater (4) completed the CME for this substance taking 
into account the mobility the titanium ion. 

The present publication gives different generalization. mainly 
characterized the assumption that the polarizabilities the particles 
depend the direction the field that, except for very special cases, the 
crystal anisotropic, i.e. its dielectric constant depends also the direction 
the field. Further assumptions are that all particles are alike, that the 
principal axes the tensors polarizability are oriented parallel each 
other, and that the particles are not too packed. These assumptions 
will specified later. 

The new equations can applied, the first place, artificial dielectrics 
used lens material for microwave frequencies. Another application can 
made crystals chemical elements. 

the following sections the term ‘‘anisotropic used when the 
polarizability the particle tensor, the term ‘‘anisotropic 
used when the dielectric constant tensor. The anisotropy the crystal 
due the anisotropy the particles the structure crystal 
lattice, general, both them. 


received original form December 11, 1952, and revised, August 
from Eaton Electronics Research Laboratory, McGill University, Montreal, 
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to 


ISOTROPIC PARTICLES CUBIC LATTICE 

giving short review how the CME usually derived for crvstals 
with cubic structure and isotropic particles lay the foundation for the 
following section anisotropic crystals. 

produced. 

The field the dielectric definition 


or 

there are particles per unit volume, then 
Np. 


Once the relation between and known, can computed from 


The field produced the point y,s, dipole parallel 


isotropic crystal the directions and are the same. nothing 
mentioned the contrary the field will considered parallel the 
rectangular coordinate system. 


ir 
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Fig. represents dielectric between the two condenser plates and 
The field produced the electric charges and The average 
field between and modified the dipoles such that 


with its center particle The radius the sphere chosen 
large compared the average distance between two neighboring dipoles. 
The distance between the plus and the minus charge one particle must 
small compared allowing the physical dipoles finite size replaced 
point dipoles and allowing moments higher order neglected. 

The field acting the particle the sum three terms: 


where the field the charges the condenser plates and 

the field the dipoles within 

the field all the other dipoles. 
assume coordinate like that Fig. with its the 
point Fig. The coordinates dipole within the sphere may 

The field all the dipoles inside the sphere 


Under the assumption that the particles are arranged cubic lattice, then 
for geometrical reasons: 


k k 
fore, 


The dipoles outside can dealt with though the space 
uniformly density dipoles per unit volume because their distances 
from are large compared the distance between two neighboring dipoles. 
One easily obtains 


From [10], [12]: 


Using [3] 
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From [5] and [14] find: 


[15] 3 = 
This the Clausius-Mosotti equation (CME) which can also written: 


[16] 
ANISOTROPIC CRYSTALS 

this section the assumption isotropy dropped. Both the polariz- 
ability and the dielectric constant are tensors. Even the field applied 
such direction that the polarization parallel the field acting 
particle may have components normal 

anisotropic dielectric there exist three mutually perpendicular 
directions such that: 


assume dielectric with only one type particles whose polarizability 
given 
br = + Ory Eaip.y + 


From [3] follows that: 
The dipoles have component the The field 
not have components these directions. 
Again the field computed drawing sphere and splitting 


can computed from [6], [8]: 
with 
the sums are extended over all 


= 


except the dipole the center 


CLAUSIUS-MOSOTTI EQUATiON 


Introducing [10], [20], [21], and [19] into [18} gives 
From and [4] 


[24] E, + 3 4nP, 1) Pz Libs, 


Substitution [24] reduces [23] to: 


From these equations can computed: 
[27] = A; = Oy 2Azy)/laix'. 
Also B,, and can computed from [26]. 


Six more equations result from [26] interchanging the and direc- 
tions. B,,, and are computed from them well, find: 


Therefore 


Using the abbreviation 


find from [25] and [28]: 


This the generalization the CME [15] and [16] for anisotropic dielectrics 
with one type particles. 

may called ‘‘apparent contains the effect the 
crystal structure represented and the effect the particles’ aniso- 

The effect the crystal structure can eliminated adding [28] the 
two other equations resulting from [28] when replaced and 

Since 
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find: 

with 


Equation [33] may also called generalization the CME [15] for aniso- 


tropic dielectrics, being kind average polarizability. 


Remark 
the main axes the tensor polarizability are parallel those the 
tensor dielectric constant, 


then 
[35] 
A 3 Arr 
the particles are isotropic, a:; that the anisotropy 
the crystal results from its structure alone, then 
136] As A = a. 


Equation can also obtained from van Santen and (6) 
when the latter specialized for one particles. 
Remark 

Another special case isotropic dielectric with anisotropic particles: 


reduces 


e+2 
reduces 


Comparing [37] with [38], with [29] and [34], gives 


[39] the condition that isotropic dielectric consists anisotropic 
particles. Since, usually, the are functions the frequency the applied 
electric [39] can fulfilled only for one special frequency. 
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If, however, the dielectric and the particles are isotropic, equation [39] 


Remark 

The meaning equations like [39] may seem uncertain because 
the sphere limiting the sums [22] has not been defined distinct 
manner. However, are mainly determined those terms for 
which small, whereas the sum the terms with larger rapidly tends 
zero. Therefore, the uncertainty can removed extending the sums 
infinity. 


APPLICATION ARTIFICIAL DIELECTRIC 


artificial dielectric may consist thin circular metal disks arranged 
and the flat particles lie parallel the plane, then 


These are the limiting values for the polarizabilities oblate ellipsoid 
vanishing thickness The particles may arranged rectangular lattice 

the constants are the same all three directions, the struc- 
tural contribution the anisotropy the crystal zero, A,=A,=0 

been computed from [22] for three different radii, r’, spheres (Fig. 2): 


0.13759 for 50. 


The fact that the three values obtained for are almost the same confirms 
that the near neighbors particle contribute predominantly the value 

From [32] and [41], 

assume particles radius and thickness 0.02. The polariz- 
abilities are computed from [40]: 


[43] 0.0267. 


From [35] 

Bz = Be 
0.0268. 


This result confirms Siisskind’s (5) statement based Cohn’s (1) experiments 
that the dielectric constant smaller than could expected from the 
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CME and that approaches the values predicted the CME 
c/a tends unity. 

When the formulas are applied the magnetic field must borne 
mind that the polarizabilities are negative. follows from measurements 
Peppiatt (3) that this case also the theory agreement with 
experiments. 
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SPECTRUM THE AIC 


ZEEMAN? 


ABSTRACT 


extensive system bands the visible and near infrared regions the 
spectrum has been observed emission when aluminum present King 
Furnace. These bands have been attributed the AIC molecule. metallic 
carbide spectrum has been described before. vibrational analysis the bands 
has been made. The 350.01, kaysers were obtained. 


INTRODUCTION 


While exposures were made the violet AIH band system with King 
Furnace source, new band system was observed the visible region 
the spectrum. The system consists number remarkably long sequences. 
The bands are degraded the red. 


EXPERIMENTAL 


(1) The King Furnace 

The King Furnace used was designed the writer and built the N.R.C. 
Laboratories, Ottawa. Some ideas the design were obtained from previous 
workers (1, 3). 

Fig. the furnace shown section. The graphite tube which served 
element was clamped its ends two split graphite cones. 
These cones fitted into two brass cones which were cooled water running 
through outer jacket, care being taken avoid the formation ‘‘dead 
the running water. The brass cones were hard soldered brass 
flanges. These flanges were bolted brass cylinder, which housed the 
insulation materials. 

The thermal insulation the furnace element consisted layer 
Zirconia Insulation Grain (or granulated charcoal) followed 
laver Johns Manville Firecrete castable refractory. These insulation 
materials formed two concentric and were bounded the inside 
graphite tube and the outside perforated steel This thick- 
walled formed unit its own and could slide into the 
main brass cylinder. 

assemble the furnace the was first put into position 
inside the main brass cylinder. The flanges carrying the brass cones were then 
bolted the main cylinder. graphite furnace element could now inserted 
through the central opening the furnace and its ends clamped into position 
means the split graphite cones. These graphite cones were held firmly 
position means two steel plates which pressed against the outside surfaces 

Manuscript received July 14, 1953. 
The work described the present paper was carried out the writer while was 
holding postdectorate fellowship the National Reseurch Council Laboratories, Ottawa, Ontario, 
Canada. The manuscript was the Merensky Institute for Physics, Stellenbosch, 
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the cones shown Fig. The graphite cones were threaded the 
outside end take pull-out device. burnt-out tube could replaced 
about fifteen minutes, provided the furnace was cold. 

The outer ends the brass cones mentioned above were closed means 
two brass plates bolted against O-rings the walls the cones. Quartz 
windows were waxed the central portions these plates and were cooled 
running water. Two metal disks could rotated front the windows 
prevent the deposit evaporating materials. This precaution was found 
necessary especially when tube was heated for the first time. 

The power for the furnace was supplied kva. step-down transformer 
with large Variac the primary circuit for control the temperature. 
The secondary terminals were connected the flanges carrying the brass 
cones. Electrical insulation between the brass flanges and the main 
was obtained means bakelite rings and washers and rubber gaskets 
shown the figure. 

The graphite furnace element was tapered shown the figure obtain 
longer region uniform temperature. The furnace was designed for absorp- 
tion work which long path advantageous. The life tube below 
2600°C. several weeks; but decreases perhaps one day 2800°C. and 
about hour 3000°C. This decrease lifetime due the rapid 
evaporation graphite the elevated temperatures. The life tube also 
depends the pressure inside the furnace and the tvpe gas with which the 
furnace filled. 


(2) Technique for Exciting Carbide Spectra 

was found that when substances like silicon, boron, and aluminum were 
placed inside the graphite furnace element their carbides were formed 
elevated temperatures. Accordingly the graphite tube either burned through 
became nonconducting large extent. either case the furnace element 
became useless. order avoid this, second and shorter graphite tube 
was inserted liner inside the tube which served the furnace element. 
The substance interest was placed inside the liner which turned into the 
corresponding carbide when the furnace was heated. this way carbide 
liner was made situ. Commercial carbide liners were found unsuitable 
because the cement which held such carbide tubes together distilled off 
elevated temperatures, with the result that the tubes crumbled away. 

While experimenting with aluminum aluminum carbide liner 
hydrogen atmosphere the violet AIH bands were observed with extreme 
intensity. Besides the AIH bands new system closely spaced bands was 
observed the visible region. Under low dispersion the bands appear single- 
headed, but the third order some band heads split into very close doublets. 
The bands are degraded the red and form very long sequences. 


(3) Evidence Favor AlC the Emitter the New System 


The following experiments were carried out order discover the carrier 
the new spectrum. new graphite liner was outgassed vacuum 
2300°C. The furnace was then filled with hydrogen pressure cm. 
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raising the temperature 2150°C. the new spectrum did not appear. 
When aluminum was placed inside the liner the hydrogen atmosphere the 
new bands appeared strongly 2150°C. With aluminum the liner, the new 
bands also appeared helium and argon atmosphere 2150°C., but only 
very weakly nitrogen atmosphere. This might due the rapid formation 
aluminum nitride the high temperatures. The aluminum used was the 
form wire. This wire was supplied for aluminizing mirror surfaces and was 
high degree purity. 

The above experiments indicate that aluminum essential for obtaining 
the bands. the next experiments carbon was introduced under controlled 
conditions. zirconium oxide liner was carefully outgassed vacuum 
high temperature. Aluminum was then inserted into the liner and the furnace 
filled with argon pressure cm. 2150°C. some bands were 
observed, but not the new spectrum. short piece specpure carbon rod 
was then placed inside the liner well fresh aluminum, and the furnace was 
again filled with argon. 2150°C. the new bands were developed strongly. 
The above experiments indicate that the new bands are only excited when 
aluminum and carbon are present the furnace liner. Hence AIC the most 
likely emitter. 

VIBRATIONAL ANALYSIS 

The bands were photographed the first, second, and third orders the 
grating spectrograph Ottawa. This concave grating has 15,000 lines 
per in. The spectrum iron was used the comparison spectrum. 
the vibrational first order exposures taken III plates were 
measured Gaertner comparator. Table the Deslandres table the 
new band svstem given. The intensities and spacings the bands are such 
that the (0,0) band could easily recognized. The long sequences also 
the construction the Deslandres table. Fig. spectrogram 
showing some the bands reproduced. 

order evaluate the vibrational constants use was made the expres- 
obtain the following relations: 


and 
{2] => (2woXo 6Gwovo) Gwo Vor. 


Hence graph against should give straight line with slope 6woyo. 

against should give straight line with intercept (wo and 
slope Such graph plotted Fig. for the lower state 
using already obtained. vields 347.98 and 2.038 
similar manner one obtains for the upper state wo’ 277.96 
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Fic. determination and for the lower state the observed AIC bands. 


Using these vibrational constants one finds that the heads the AIC bands 
can represented 
order obtain the dissociation energy the lower state, was 
plotted against and extrapolated. plot against was also made. 
The results obtained are rather uncertain, but estimated 1.6 0.2 ev. 


TABLE 


VIBRATIONAL CONSTANTS THE MOLECULE THE UPPER AND LOWER STATES THE VISIBLE 
BAND SYSTEM 


Constant Lower state Upper state 
von | 17234.0; kaysers 
2.022 kaysers 0.831 kaysers 
2.038 kaysers 0.846 kaysers 
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ROTATIONAL ANALYSIS 


The (0,0), (0,1), and (1,0) bands were photographed the third order 
high contrast plates. These are the only bands which are not badly masked 
the overlapping adjacent bands. Each band has two heads about 
equal intensity, very close together. Groups lines regularly spaced seem 
form series, but the lines each group are not sufficiently resolved establish 
the number branches each band. The second difference these groups 
lines very small, about 0.033 kayser. Thus far complete rotational 
analysis could not carried out and therefore not possible ascertain 
what type transition the band system represents. 


CONCLUSION 


view the fact that the bands have been observed thermal emission, 
the two electronic states involved the new bands must low-lying states 
and appears quite probable that the lower state the ground state the 
AIC molecule. the other hand, difficult account for such small 
value the vibrational frequency the ground state when that 
AlO AIF their respective ground states about twice large. For 
unambiguous assignment AIC, investigation the isotope effect would 
desirable. However, the experimental evidence for AIC presented above 
seems fairly conclusive. Work under still higher resolution will required 
order satisfactory rotational analysis this band system and 
determine the nature the upper and lower states. 


ACKNOWLEDGMENTS 


The writer wishes thank the National Research Council Canada for 
the fellowship awarded him. The writer also wishes thank Dr. Herzberg 
who suggested the problem, for his guidance and valuable suggestions the 
course the work. 


REFERENCES 


HERZBERG, Molecular spectra and molecular structure. Vol. Spectra diatomic 
molecules. Van Nostrand Company, Inc., New York. 1950. 96. 

Trans. Am. Electrochem. Soc. 56:97. 1929. 


i 
é 
3 


BIREFRINGENCE CRYSTALS AND THE 


ABSTRACT 


Propagation plane electromagnetic waves homogeneous ionized gas 
uniform magnetic field compared with the propagation light optically 
inactive birefringent crystal. well known that propagation crystal may 
described using system real orthogonal axes for which the dielectric 
constant given diagonal matrix. This paper shows that propagation 
plane waves the ionosphere may described similar manner, the medium 
having effective dielectric constant given diagonal matrix, provided that 
system orthogonal axes used for the description the compo- 
nents the field vectors. This set component axes (which quite different from 
and not confused with coordinate axes) equivalent resolving the field 
vectors into components parallel the magnetic field and two contrarotating 
circular components plane perpendicular the magnetic field. expression 
giving the velocity each the two modes propagation given direction 
and expressions for the amplitude each component the field vectors are 
obtained (equations and 44). Provided that one accepts the concept 
complex velocity propagation, the results hold when electron collisions are 
included. When electron collisions are neglected, possible form double- 
sheeted surface, called the normal velocity surface, which some assistance 
visualizing the manner which the velocity propagation the plane waves 
each mode changes with direction. 


INTRODUCTION 


comparison double refraction crystals and the ionosphere has been 
made Lange-Hesse (6). Using the theory crystal optics showed that 
plane wave optically inactive anisotropic crystal was characterized 
two velocities propagation, thus explaining the birefringence the medium. 
then obtained corresponding equations for the propagation plane waves 
collision-free homogeneous ionized gas uniform magnetic field and 
compared them with those obtained for the crystal. similarity existed since 
principal axes coordinates could found for both media (relative which the 
effective dielectric constant could written its simplest form diagonal 
matrix). also existed since both media were birefringent—but the 
details the were very dissimilar. For example, might 
thought that the ionized gas the magnetic field could regarded medium 
with physical properties which are symmetrical about axis parallel to. the 
magnetic field, that the gas might expected have properties similar 
those shown crystals with electrical axis symmetry (uniaxial crystals). 
actual fact the properties are quite different, since firstly, each mode 
propagation linearly polarized crystal, but elliptically polarized the 
ionosphere; and secondly along the axis symmetry the velocities the two 
modes propagation are equal (uniaxial) crystal but are not equal the 
ionosphere. 

believed that propagation the ionosphere could more closely com- 
pared with propagation optically active anisotropic media, than with optically 

March 31, 1958. 
from the Defence Research Telecommunications Establishment, Radio Physics 
Ottawa, Canada. Work performed under Defence Research Board number 
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inactive anisotropic media, since elliptical polarization always associated 
with each the two modes propagation the latter. However the interpre- 
tation the physical characteristics propagation optically inactive 
anisotropic medium simpler than that for optically active anisotropic 
medium. the present investigation therefore, like Lange-Hesse will still 
carry out comparison propagation the ionosphere with that 
inactive crystal, but introduce the added complication damping due 
collisions electrons with other particles. found that the ionosphere 
(regarded slowly varying medium) the propagation plane wave 
always characterized three components the field vectors, one which 
linearly polarized along the direction the uniform magnetic field and the other 
two circularly polarized the planes perpendicular the magnetic field. 

Throughout this paper the rationalized m.k.s. system units, the time factor 
and the representation plane waves the factor 
are used. The derivation the effective dielectric constant 


PROPAGATION ANISOTROPIC MEDIA 


well known that the field vectors representing electromagnetic field 
satisfy the equations (Ref. 14, 23) 


[3: 

crystal 

B = 

and 


where are the permeability and permittivity free space and are indepen- 
dent the directions (D, respectively, and e,, the relative 


permittivity the crystal, tensor. 
the ionosphere, which will considered ionized gas uniform 


magnetic field, 


[8] B = 
[9] 


where the oniy new constant the conductivity the medium, tensor. 

From examination the set equations and the set 10, there 
would appear connection between the propagation electromagnetic 
waves the two media. But introducing the effective displacement vector, 
and writing 


Ey 
§ 
4 
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equation some similarity apparent, since for both media one now has 


and since divcurl one also has 


The equations 12, 13, are, course, true for any medium. Next since the 
field vectors have time factor one obtains from 12, 13, 


Then using one observes that crystal 

and the ionosphere 


where and are tensors. These constitutive equations and 
are the kev the properties the two media.* 

restrict the mode propagation plane waves, then each the 
vector field components will contain the factor 
where the vector from the origin point the field, and isa unit vector 
perpendicular the planes constant phase, that with respect any set 
orthogonal axes Oxyz, the direction cosines are (a, The phase 
velocity and will seen shortly dependent the direction 
Using 20, equations 15, 16, become 


[23] 


and eliminating 


9 = 5 = 


*Some authors prefer the use the polarization vector equation 19, which would then 
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Some important features the propagation are seen immediately from these 
equations (Fig. 1), for evident that and form right orthogonal 
system, and that the plane and (Ref. 372). 


Instantaneous position the effective displacement, electric, and magnetic field 
vectors and anisotropic media with respect direction propagation (Note 
that and form orthogonal set.) 


anisotropic crystal identical the displacement vector and the ionosphere 
effective displacement vector. 


Next, equations 21, 22, and will used discover more about the velocity 
propagation and the relationships between the components the field 
conducting anisotropic crystal will studied first, and then comparative 
study made propagation ionized gas uniform magnetic field. 


PROPAGATION NONCONDUCTING ANISOTROPIC CRYSTALS 


well known that nonconducting anisotropic crystal possible 
orient orthogonal system axes that and are connected 
diagonal matrix containing real elements, i.e. 


The axes Ox, Oy, are then called the principal electric axes the crystal. 
Using 24, and 18, one obtains 


and resolving this into its vector components (using 25) one finds 


a(s. 


[27] 
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COMPONENT 


COMPONENT 


(a) 


(b) 


Fic. features propagation ionosphere. 
(a) Principal components electromagnetic field vectors. 


The normal velocity surfaces, and the polarization the effective 
‘or propagation plane waves along cardinal axes Ox, Oy, and Oz. 


(c) Polarization the electric vector for propagation plane waves aiung Ox, Oy, and 


The two modes propagation are propagated independently, and the 
existence one does not require the existence the other. (Coupling between 
the two modes does exist, however, when the ionized gas not homogeneous, 
Ref. 11.) 

For the hollow cone which includes all directions propagation inclined 
constant angle the direction the magnetic field, the velocities the 
iwo modes propagation are constant (equation 43), and hence this cone 
‘equation 44) the (¢) components the field vectors parallel the magnetic 
remain unchanged amplitude and phase, and the two circularly polarized 
and components the field vector change phase but not amplitude. 
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Normal Velocity Surfaces 


normal velocity surface surface such that the radius vector from 
the origin everywhere proportional the phase velocity the plane waves 
propagated the direction the radius vector (see section 3). 

Corresponding the normal velocity surface for crystal, shown Fig. 
there normal velocity surface for propagation ionized gas uniform 
magnetic field. for double sheet, one octant which shown 

These figures give true physical representation phase velocity any 
direction when the collision frequency negligible, for, using the relationships 
given appendix 30, 31, 32, 33, may shown that when the frequency 
the propagated wave (w) and the critical frequency satisfy 


Since the velocity independent the strength the uniform magnetic 
field, the sheet PQR corresponds the so-called ‘ordinary’ mode propagation, 
and the sheet UVW the extraordinary mode (Ref. 10). The arrows attached 
to, and around the axes Figs. 4(b) and 4(c) give the polarization the 
effective displacement vector and the electric vector respectively for electro- 
magnetic plane waves propagated the direction the axes. 

These directions polarization the effective displacement vector are 
obtained from equation 44a after using equation express terms 
D,, similarly, the directions polarization the electric vector are 
obtained from equation will seen that, for propagation along the mag- 
netic field, the ordinary mode possesses only the circular component, and 
the extraordinary mode possesses only the circular component. For 
propagation perpendicular the magnetic field the ordinary mode 
possesses only the linear component. However the extraordinary mode 
possesses only the two circular components which, being equal amplitude 
interesting recall Fig. which was pointed out that general, has 
component the direction propagation. may (Fig. that 
this component the direction propagation vanishes for both modes 
propagation along the uniform magnetic field, and for the ordinary mode 
propagation perpendicular the uniforni magnetic field, but does exist for the 
extraordinary mode propagation perpendicular the uniform magnetic field. 
The polarization the current vector similar that illustrated Fig. 
for 

From the expressions 30-33 appendix which are dependent the value 
(and therefore the electron density) will seen that the 
shape the normal velocity double sheet must change wave propagated 
through the ionosphere. The nature this change shape the sheet may 
followed with the help Figs. and (similar the diagrams Refer- 
ences and 11) which show the variation (c/v)? (or square the refractive 
index) with for both modes and for any direction propagation but neglect 
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Fic. Variation (refractive with (critical for plane waves whose 
direction propagation inclined any angle the direction the magnetic field. (Note 
that collision frequency has been neglected and that wo= eB/m negative for electrons.) 
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collisions. For example, when the frequency propagation greater than the 
gives the characteristics the propagation any direction 
point the ionosphere having the critical frequency given the intercept 
with the abscissa. The velocities corresponding the strips 
and Fig. give the velocity variations the ordinary and extra- 
ordinary modes (OR OQ) and OV, Figs. the direction 
propagation changes from along perpendicular the direction the uniform 
magnetic field. one ascends through the ionosphere, the critical frequency 
initially increasing, the line moves the right, and the normal 
velocity double sheets change their shape accordingly. For instance the line 
plane through the direction the magnetic field shown Fig. will 


ORDINARY 


MODE 


ATTENUATION 


Section the normal velocity surface point the ionosphere where the 
critical frequency corresponds (Fig. 5(a)). Note attenuation corresponding 
negative values c?/v*. 


observed that either mode may the inner outer sheet depending the 
critical frequency (i.e. electron density) and also that when (c/v)? imaginary 
(as the attenuation which represents may also illustrated (as 
the shaded section). 


CONCLUSIONS 


The propagation plane electromagnetic waves homogeneous ionized 
gas uniform magnetic field has been compared with that optically inactive 
nonabsorbing anisotropic crystals. 

has been confirmed that normal velocity surfaces (consisting two sheets, 
one corresponding each mode propagation) may drawn for both media. 
optically inactive crystals well known that these surfaces generally meet 
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four singular points, defining two optical axes along which plane waves 
associated with each mode are propagated with equal phase velocities. has 
been shown that the ionosphere, when collisions are neglected, these surfaces 
not meet any point, and, given direction, the velocities propagation 
associated with each mode are general different. 

has also been established that for propagation the ionosphere even when 
electron collisions are included, the electromagnetic field vectors are very 
simply related each other diagonal matrix transformations, the field 
vectors are resolved into three special components, the first which linearly 
polarized along the direction the earth’s magnetic field, and the second and 
third which are two contrarotating circularly polarized components lying 
the plane perpendicular the magnetic field. Here differ from Lange-Hesse 
(6) who reported that for collision-free ionosphere the second 
components were polarized. also differ from Saha and Banerjea 
(Ref. 12, equation 22) who, without reference propagation crystals, had 
suggested three components similar ours, but implied that they travelled 
with three velocities. The latter authors made reference the diag- 
onal matrix property. 

Apart from the simplicity the diagonal matrix transformation connecting 
the field vectors when they are resolved the special manner suggested above, 
one also obtains simple picture the polarization the field vectors when 
the plane electromagnetic wave propagated along perpendicular the 
magnetic field. From the general expressions equation 44, the directions 
polarization (say) the E-vector shown Fig. for propagation the 
direction any the axes and for either mode, are simply obtained. 

the plane wave passes through the ionosphere, the electron density, and 
hence the critical frequency w,, changes. Provided the ionosphere may still 


regarded slowly medium, the concept normal velocity surfaces 


still holds but these surfaces change their shape conformity with the charac- 
teristics obtained previous authors and illustrated Fig. 


APPENDIX 


The Conductivity and Effective Dielectric Constant Matrix for Ionized Gas 
Uniform Field 


This matrix has been obtained previously several authors (Ref. and 
Ref. 12, equations and 12). Its derivation repeated here, following Hines 
(Ref. 4), but using the m.k.s. nomenclature order provide convenient 
reference for the form the effective dielectric constant matrix 


= BO 


given equation 30. 
The equation motion charged particle mass and charge 
constant magnetic field 
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where velocity, 
and electric and magnetic field (due the electromagnetic 
radiation), 


collision frequency, 
permeability free space. 
now assume that the field vectors associated with the electromagnetic 


radiation have common factor 


that the radiation plane wave with phase velocity and direction 
propagation defined direction cosines (a, the right orthogonal system 
axes then the term containing usually neglected since 
generally written 


usual consider only the motion electrons, since the velocity other 
(heavier) charged particles much smaller. Hence the current NeV, 
where the electron density, related the electric field the expression* 


Now 
permittivity free space, 
2 Ne’ 
Re 
(jw+ 


where are the components parallel the coordinate axes 
then one obtains from 


the movement electric charges considered terms displacement rather than velocity, 
one could appropriately use the polarization vector where 
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this, also evident from equation A-3, the effect the collisions 
changes the conductivity tensor from v). Then inverting equa- 
tion A-7, the conductivity tensor given 


2 2 
2 2 
{.\-8] = W1W2 — a + We WoW3 + E. 
W1W3 + awe @oW3 — AW) a 4 


now, special coordinate axes O(x, are chosen that the axis lies along 
that 


Hence the displacement vector and the electric field vector 
related equation are connected the expression 


2 2 
2 
jwa 


which may written the form quoted equation 30, namely 
APPENDIX 
Reference List Symbols, etc. 


magnetic, induction, and current density field vectors; each contains the 
plane wave factor 


p—electric charge density. 

propagation the plane wave. 

velocity and angular frequency the plane wave. 


right orthogonal real axes. the crystal these are 
the principal electric axes, and the ionized parallel the uniform 
magnetic field. 
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21. 
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system right orthogonal complex component axes (com- 


plex principal electric axes) the ionized gas the uniform magnetic field. 
permittivity free space. 
constant, permittivity relative free space crystal 
and ionosphere respectively; these are matrices. 


1 2 


Ou; 
the coordinate axes. 


E,, components the field vectors etc. resolved 
the directions the real axes Ox, Oy, Oz. 


resolved the directions the complex ‘component’ axes On, Of. 
associated with principal dielectric constants crystal. 


velocities the two modes propagated the direction 


(a, 
propagation. 


associated with principal dielectric constants 
ionized gas uniform magnetic field, given etc. 

When there are collisions these components are real imaginary, and 
when there are collisions they are complex. 


2 
aw, 


e—mass and charge electron. 

N—electron density ionized gas. 
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COINCIDENCE STUDIES THE DECAY INDIUM 


AND 


ABSTRACT 


coincidence circuit resolving time sec. has been used examine 
the decay Coincidence experiments are described show that the 
five gamma-rays which follow the K-capture process the 0.722 and 1.300 Mev. 
radiations are both cascade with the 0.556 Mev. gamma-ray and that the spins 
and parities the excited levels energy 0.556, 1.278, and 1.856 Mev. 
respectively. The observed coincidences 
may ascribed combination inner bremsstrahlung associated with the 
strong Mev. transition and weak group intensity less than 
beta per disintegration and end point about 0.8 Mev. The positron spectrum 
has branching ratio end point between 1.0 and 1.4 Mev. 


INTRODUCTION 


The 50-day activity has been studied many workers. The main 
features the strong spectrum were established Lawson and Cork (20) 
and confirmed others. The 50-day activity due isomeric transition 
from excited state energy 0.192 Mev. and spin the ground 
state spin (6). The spectrum due allowed transition from the 

addition this dominant mode decay, K-capture leading excited 
state was discovered Maienschein and Meem (24) and their results 
confirmed and extended Boehm and Preiswerk (4) and Mei, Mitchell, and 
Zaffarano (25). This process, intensity the intensity, leads 
excited state energy 1.28 Mev. and de-excitation occurs through 
the emission cascade gamma-rays energies 0.722 and 0.556 Mev. The 
presence weak radiation energy around 1.3 Mev. has also been established 
and its role crossover transition suggested. Boehm and Preiswerk (4) de- 
tected weak component the decay process and estimated its branching 
ratio 0.01% and its end point 0.650 Mev. 

This work was originally undertaken determine the spins the excited 
states angular correlation measurements. The difficulty 
terpreting the results the angular correlation measurements (14, 15) and the 
equivalent measurements Steffen (28) led re-examination the decay 
scheme. The discovery Steffen and Zobel (29) that the 1.278 Mev. level 
fed from the ground state has been great assistance in- 
terpreting the results. Precise measurements the gamma-ray energies carried 
out this laboratory (17) showed that the 1.300 Mev. radiation could not 
the crossover transition originally proposed. addition, these measurements 
indicated that this radiation was cascade with the 0.556 Mev. radiation. 
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The coincidence experiments described this paper support the decay 
scheme Fig. already proposed the basis energy measurements (16), 
set upper limit the number negatrons proceeding excited state 
assign spins the excited states and determine the yield and 
approximate end point the positron component. 


SN ii¢ 


ALL ENERGIES 
MEV. 


1.278 


L271 


Fic. decay scheme for 


EXPERIMENTAL METHOD AND RESULTS 


The coincidence measurements were carried out with circuit resolving 
time three five millimicroseconds similar design that described 
Bell, Graham, and Petch (2) which the reader directed for detailed 
account its performance. The gamma-rays were detected with anthracene 
were strapped either 931A photomultiplier tubes. The entire tube and 
crystal was covered with aluminum foil shield which was connected the 
cathode the photomultiplier. Further shielding Scotch electrical tape was 
provided ensure that the counters were light-tight. For gamma-gamma 
coincidence measurements the crystals were provided with lead shields mm. 
thick eliminate betas and attenuate the strong 0.192 Mev. gamma radiation 
from For beta counting the crystals were covered with mgm./cm.? 
aluminum foil. Each photomultiplier was mounted small portable copper 
box containing its resistor network and preamplifier. The photomultipliers were 
operated 2000 

The performance the circuit was tested obtaining resolution curve 
under the conditions operation the experiment. With clipping cable 
cm. long, the circuit gave efficient detection with resolution times 
sec. for anthracene and sec. for stilbene good agreement with the 
experience Bell, Graham, and Petch. (Since the circuit possessed resolution 
time 0.8 sec. when operated from square wave pulse generator (10) 


\96.2% 
556 
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whose rise time was 0.5 sec., clear that the largest contribution 
the resolving time the circuit found the properties the phosphor 
used.) The resolution curve was obtained measuring the coincidence rate 
function the delay introduced between one the scintillation counters and 
the fast coincidence detector. The following simplified version the delay line 
described Bell a/. was used introduce this delay. The fast pulses from 
the two counters were fed through 100 ohm cable (RG/7U) the two ends 
100 ohm line 2.4 meters long formed with solid central copper conductor and 
copper shield U-shaped cross section. The conductor was held place with 
small pieces lucite. The fast coincidence detector and its associated pre- 
amplifier Los Alamos model 500) rode this line and the pulses entered the 
coincidence unit through phosphor bronze spring which made good contact 
with the copper conductor. Since shift the contact point increased the 
length the cable leading one counter and reduced for the other, change 
contact position cm. corresponded delay sec. Thus this 
variable delay line allowed total delay 1.6 sec. introduced. 
Larger delays were provided changing the length the cable from one 
the counters the variable delay line. The velocity pulses cable 
given the manufacturer 0.83 while the velocity the variable line 
the velocity light. 

After checking the resolution, the contact was set the point corresponding 
the peak the resolution curve and the experiment begun. Each measure- 
ment involved determination the difference between the coincidences 
short interval time, min., and the number chance coincidences 
equal interval. The chance coincidences were obtained switching meter 
length cable one other the leads from the counter the coincidence 
detector. The chance rate was corrected for the small effects attenuation 
the delay cable. The proper correction factor was found exciting each photo- 
multiplier with separate gamma source and examining the change the 
chance coincidence rate with the length the delay line used. For the meters 
cable used for measuring the chance rate the correction factor for attenuation 
was 1.15. 

general, was not necessary cool the photomultiplier. Although the 
noise contribution was very high the high voltages used, the resolving time 
the coincidence circuit was sufficiently short that the noise coincidence rate 
was completely negligible. In-the study the coincidence rate 
function lead absorber thickness, one crystal was exposed the full 
intensity the source while the other was shielded with large amount lead. 
this case, was found that coincidences between noise pulses from the shielded 
counter and the very large number pulses from the unshielded counter were 
sufficiently frequent troublesome. For such experiments, the shielded 
photomultiplier and its crystal were cooled with liquid air with consequent 
reduction the noise coincidence rate factor 10. 

Coincidences between the 0.556 and 1.300 Mev. Gamma-rays 

From the gamma energy measuremencs already reported (17), appeared 

probable that the 0.556 Mev. radiation was coincidence with both the 0.722 
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and 1.300 Mev. radiation shown Fig. The following coincidence ab- 
sorption experiment was carried out test this part the decay scheme. 

Anthracene scintillation counters shielded with mm. lead were placed 
opposite sides the source distances and cm., respectively, from the 
source. Lead was placed between the source and the distant counter and the 
coincidence rate measured function the thickness lead. The coincidence 
rate with source 100 mc. was about 200 c.p.m. with absorber 
present and dropped factor roughly 1000 with cm. absorber 
place. The signal chance ratio for zero absorber was six. With the shielded 
crystal and the photomultiplier immersed liquid air, the signal chance ratio 
was still two with cm. lead screening the cooled counter. 

order minimize the effect drifts the electronic circuits, counts were 
taken alternately with and without the absorber place and the measurements 
for absorber normalized the value 1000. 

Since the measurements were not made under conditions geometry, 
was impossible identify the energies the gamma-ray components from 
the absorption curve reference the known absorption coefficients. 
method identification, the experiment was repeated with source 
precisely the same geometric arrangement. 

The experimental results are shown Fig. Curves and represent the 
data while curve represents the data. Curve was taken with the 
side-channel discriminator the shielded counter set low value while 
curve was taken with this discriminator set discriminate strongly against 
the 0.556 Mev. radiation. All data have been corrected for background coinci- 
dence rate 0.3 the normalized scale. apparent that both curves show 
two distinct components, the softer which, common both, corresponds quite 
well with the theoretical curve for equal mixture 0.556 and 0.722 Mev. 
radiation and the harder which seems energy slightly greater than 
that the Co® radiation mean energy 1.2 Mev. 

The decay accompanied six gamma-rays (16, 17) whose energies 
and relative intensities (given brackets) follow: 0.192 Mev. (5.4), Mev. 
(1), 0.722 Mev. (1), 0.576 Mev. (0.03), 1.300 Mev. (0.04), and 1.271 Mev. 
(0.009). interpretation the absorption curve must based these 
energies and intensities. The 0.192 Mev. radiation due isomeric transition 
the 72-sec. ground state that cannot coincidence with the 
other listed above. order decrease the chance coincidence rate, 
the effective relative intensity this radiation was reduced less than 0.2 
means lead shielding and the discriminator settings the coincidence circuit. 

The strong component the coincidence absorption curve has slope con- 
sistent with equal mixture 0.556 and 0.722 Mev. radiation and accounts 
for almost all the gamma-gamma coincidences observed. must therefore 
identified due the Mev. cascade observed many workers 
(4, 24, 25). 

The weak component must ascribed coincidences between either both 
the 1.271 and 1.300 Mev. radiations and one more the 0.556, 0.576, and 
0.722 Mev. gamma-rays. support this assumption, discrimination against 


i 
| 
eGOo0uoU 


JOHNS ET AL.: DECAY OF INDIUM 114 


“VE 


a 


uci 


\ 

< \ 

\ 


LEAD THICKKESS (CM) 


the 0.556 and 0.722 Mev. gamma-rays the shielded counter led large 
increase the relative intensity the high-energy component (curve B). 

Since the relative detection efficiencies the circuit function gamma 
energy are known for the conditions under which curve was taken, possible 
from the known gamma intensities compute the relative intensity the 
weak the strong component the absorption curve for all possible cascades. 
For the Mev. cascades the relative intensity 
expected about 3.5%; for the 1.271-0.722 Mev. cascades 
should about 0.5%; while for all other possibilities should less than 0.1%. 
The observed intensity under these conditions was 2.5% good agreement 
with the calculated intensity for Mev. cascade and 
quite inconsistent with any the other possibilities. 

The results these experiments confirm the fact reported several workers 
(4, 25) that the 0.556 and 0.722 Mev. radiations are cascade and addition 
show that the 1.300 Mev. radiation cascade with either the 0.556 0.722 
radiations. The existence level 0.556 Mev. established 
neutron-capture experiments (18, 30) and the presence weak 
energy 0.576 and 1.271 Mev. the decay spectrum (17) rule out the latter 
possibility. 
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Coincidences 

Although the chief mode decay involves direct transition from 
the ground state that has recently been pointed out 
Nussbaum and van Lieshout (26) that transitions the first excited state 
might expected occur small fraction the disintegrations, and 
involve spectrum with end point about 0.8 Mev., followed gamma- 
ray energy about 1.2 Mev. Their measurement sets upper limit 
the intensity this transition. The following section describes comparison 
the coincidences obtained with sources and 
attempt set upper limit the intensity such transition. 
was chosen comparison source because each disintegration consists 
(end point 0.3 Mev.) followed two cascade gamma quanta energy about 
1.2 Mev. was chosen because possesses simple spectrum end point 
(2.3 Mev.) slightly greater than Mev.) and with known gamma-rays. 

The coincidence circuit already described was used with two stilbene crystals 

detect only gamma-rays and the other with 4.0 mgm./cm.? aluminum window 
detect betas. The shielding the gamma crystal had negligible effect 
the high-energy gamma-rays interest here but reduced the contribution 
the 0.192 Mev. gamma indium and the soft bremsstrahlung radiation 
accompanying both the and decay. The beta crystal was covered with 
1900 aluminum determine its gamma counting rate. Each 
crystal subtended solid angle 0.02 the total sphere the source position. 
each the three sources the single counting rates, the rate, 
and the rate were studied. The aluminum absorption curves for 
and presented Fig. are very similar, indicating beta end point 
the region 0.8 1.0 Mev. Owing the low coincidence rates and the presence 
these absorption curves for more than 1.5 decades. 

Following the analysis described Dunworth (9) and others, 
possible compute the intensity the low-energy spectrum three 
partially independent ways from comparison the single and coincidence 
rates recorded for and The values obtained the three methods 
agree quite closely and give mean for the number low- 
energy negatrons per disintegration 

Alternatively, the observed coincidence rate may interpret- 
due bremsstrahlung. one assumes that the 
dences observed are due weak group energy about Mev. 
followed 1.2 Mev. gamma-ray, comparison the ratio for 
with that for leads intensity for this group the 
other hand, group this energy would have given ratio 
times larger than the ratio actually observed. This suggests that the observed 
beta coincidences. The magnitudes expected for such effect were calcu- 
lated for the and using the method outlined for (31). 
The number quanta per disintegration arising from inner bremsstrahlung 
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Fic. coincidence absorption aluminum for and The 
contributions have been subtracted from the and coincidences. The 
curves have been arbitrarily normalized for clarity presentation. 


with energy sufficient reach stilbene crystal covered with lead 
number sufficient account for all the coincidences observed 
both sources. The lead shielding between the two crystals was 
sufficient eliminate all bremsstrahlung coincidences arising from retardation 
the betas the beta crystal itself. 

the light this evidence, would appear that the intensity any low- 
energy group existing must less than the high-energy 
group. 

Spins the Excited States 

(i) Introduction 


Since Hamilton (13) showed that there was expected correlation between 
the directions emission the first and second gamma-rays gamma cascade, 
the experimental determination correlation functions has been used Brady 
and Deutsch (5) and other workers powerful method assigning spins 
excited states nuclei. Hamilton shows that the expected correlation function 
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pressions for the coefficients the case pure dipole and quadrupole transitions. 
The theory has been extended Lloyd (22) cover radiations any multipole 
order. The series terminates with the term unless both transitions are 
multipole order greater than two. has been pointed out Ling and Falkoff 
(21) that the competition electric and magnetic radiation transition 
may introduce interference effects which radically alter the correlation function 
expected the pure multipole theory and they present explicit expressions for 
calculating the correlation functions for mixtures. Lloyd (22) 
has shown that the interference terms introduced Ling and have 
physical reality only when the phase shift radians. the basis 
theory (3) appears that transitions mixed character can only 
expected for competition and radiations, and that even here the 
component should much weaker than the radiation. Expressions for 
correlation functions for triple cascades which the first and third 
are detected have been given Arfken, Biedenharn, and Rose (1). would 
expected, these correlation functions are much less anisotropic general 
than for two gamma cascades. 


(ii) 


The following section reports the results two separate determinations 
the correlation function for The two sources used had 
activities 100 mc. and consisted metallic indium irradiated specific 
activity 0.5 and respectively. Since the shielded 
detectors were sensitive only the gamma-rays energy greater than 200 kev., 
the effective source strength was about mc. each experiment. Despite this 
activity, the signal chance ratio was better than four one. The solid angles 
subtended the scintillation detectors the source were 0.027 and 0.045 
steradians the two experiments. Slight differences the method source 
mounting and the precautions used avoid Compton scattering existed but, 
since the individual points the two runs agreed within the expected statistical 
errors, the weighted mean the two sets data was used. 

Each source was mounted the axis laboratory prism spectrometer 
table. The anthracene stilbene crystals with their associated 1P21 photo- 
multipliers were mounted the fixed and movable arms the spectrometer. 
Each source was centered carefully observing the counting rate the movable 
detector function its angular setting. The crystal faces were covered with 
reduce the intensity the strong 0.192 Mev. radiation and 
scattered radiation, while the sides the crystals were covered with mm. 
lead for the same purpose. check the suitability the geometric arrange- 
ment, the correlation function for was determined and found agree with 
that given Brady and Deutsch (5) and others. Counting was carried 
the position, the 90° position, and chance turn until least 10,000 true 
coincidences had been recorded for each selected value 

the correlation function for the selected value small correction 
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which was never greater than 0.007 was made each value for the the 
finite solid angles subtended the crystals (7). The 180° position was also cor- 
rected for the presence annihilation radiation from the weak positron component 
present the decay. The effect this correction was lower the point 180° 
from 0.203 0.012 0.177 the increase the standard deviation 
being due the uncertainty the positron contribution. The experimental 
data, corrected described above, are presented Fig. Since the upper level 


100 120 130 140 170 180 
DEGREES 


Fic. The angular correlation function for the cascade gamma-rays. 
Solid line—the least squares fit the experimental data and the theoretical 
correlation function for the Mev. cascade modified three weak cascades. 
Dotted line—the theoretical curve for 2-2-0 cascade with mixture 96.0%, 
and zcro phase shift the first transition. 


the 0.722-0.556 Mev. cascade fed from the ground state 
spin (29), its spin not greater than while the lower level the cascade 
the ground state nucleus with presumably zero spin. Hence 
difficult conceive many cascades which both transitions could 
octupole higher order. The few which are possible give much greater aniso- 
tropy than that shown these data. was therefore felt reasonable 
fit the data Hamilton's equation, retaining only terms the points 
are given weights inversely proportional the squares their relative standard 
deviations, least squares fit the data vields the parameters 0.0028 
and 0.190. The curve described these parameters shown solid 
line Fig. Although these parameters are somewhat different from those 
obtained from the correlation measurements Steffen (28) and 
(19), the resulting curves are very closely the same. 

impossible find any pure multipole cascade which will give good 
ment with any these data. The best choice the 4-2-0 cascade which ruled 
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out the limitation the spin the upper level. Even when the possibility 
interference considered, the fit provided any the published 
sets parameters not good. The best choice for the data presented here 
with phase shift and the second pure quadrupole. This function shown 
the dotted curve Fig. very sensitive the mixture and 
radiation present. 3.5% component fits well angles near 180° but 
considerably worse angles below 150° while 4.5% component fits well 
small angles but gives too low correlation value for large values When 
considered that there are three other cascades present the decay process 
with intensities the main Mev. cascade, this measure 
agreement all that can expected. 

can assumed that the discrepancy between the experimental and 
theoretical curves Fig. represents the distortion the Mev. 
cascade the four weaker cascades, possible use the difference make 
estimate the spin the 1.856 Mev. level. The four weak cascades and 
their intensities relative the main cascade are (a) the 1.300-0.556 Mev. 
cascade (4%), (6) the triple cascade which the first and 
third gammas are detected (3%), (c) the cascade (3%), and (d) the 
cascade (0.01%). Since the main shape the correlation function 
determined the strong cascade, the spins the first three levels 
are The 1.856 Mev. level may conceivably have spin (if 
allow the possibility the upper level being fed from the 0.192 Mev. level 
relatively simple matter consider the effect the weak 
cascades, using the results the theoretical papers referred the introduction 
this section. the four possible choices, only the choice gives any improve- 
ment the agreement between the theoretical and experimental correlation 
curves. For this choice possible fit the experimental curve exactly when 
the 0.722 Mev. transition described 96.2% and 3.8% E2, with phase 
shift. Qualitatively, possible see this considering that only this case 
there strong correlation function any the weaker cascades which could 
introduce the measure distortion required. 

The fact that possible get such agreement does not constitute very 
strong proof that the spin the 1.856 Mev. level zero. Supporting evidence 
exists the absence the gamma-ray corresponding the 1.856 Mev. tran- 
sition the ground state This radiation, exists, has 
Jess than quanta per disintegration (17), value much less than 
would expected were the spin the 1.856 Mev. level 


Positron Spectrum 

The existence very weak positron component the decay was 
first suggested Boehm and Preiswerk (4). Since their measurements were 
very tentative, was felt worth while determine the and end point 
coincidence methods. Since the two 0.511 Mev. quanta accompanying positron 
annihilation are emitted opposite directions, gamma-gamma coincidences 
due annihilation radiation can observed only when the angle formed 
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the counters the source 180° and the positron yield can determined 
directly comparing the number gamma-gamma coincidences when 
not 180° with the number annihilation quanta eoincidences when 180°. 

Let a(= 0.035) the K-capture branching ratio for the cascade gamma rays 
(17) and the positron yield. Then, and represent the intrinsic efficien- 
cies the anthracene crystals for 0.511, 0.556, and 0.722 Mev. radiation, can 
1.14) correction factor made necessary because the cascade gammas are 
not emitted isotropically, while the solid angle subtended either counter 
the source (in symmetric arrangement). the ratio 
the number annihilation coincidences the number 
coincidences the 180° position. 

The apparatus used for this experiment was the same that employed for 
the previous coincidence work. order make large 
possible consistent with reasonable coincidence rates, was kept small placing 
the two crystals 23.6 cm. apart, with the source midway between them. The 
source was covered with sufficient brass stop Mev. positrons and the coinci- 
dence rate the 180° position compared with that for angle 167° between 
the counters. The rate 167° was corrected the 180° position 
using the angular correlation function described above. The effect positrons 
the correlation experiment was much less marked since the counters were 
closer the source and the source was uncovered. The ratio was 1.56 
from the known counter efficiencies, 0.00186, and 0.36 
0.03. Using these data the positron branching ratio The 
standard deviation expressed this result has been arbitrarily increased take 
into account uncertainties the branching ratio the K-capture branch and 
the counter efficiencies and solid angles. (The value previously given 
Reference (15) should have read while the value Reference (16) 
was computed incorrectly using instead 
and hence roughly four times too large.) 

principle should possible obtain measure the positron end point 
using thin source and observing the increase the number coincidences 
due positron annihilation the source covered with successive layers 
absorber. essence this the reverse the common aluminum absorption 
method determining beta end point. Since the positron spectrum was too 
weak study the spectrometer available, attempt has been made 
estimate its energy this method. 

The source indium metal was rolled into foil thickness 
and sandwiched between two mica shields thick. The mova. 
detector was swung from the 180° position the 155° position and the flat 
source oriented that this detector looked its surface from the same angle 
both positions. This symmetry was necessary prevent differences gamma 
absorption from distorting the measurements. The aluminum absorbers used 
varied thickness from 600 mgm./cm.? though the positrons were all 
essentially absorbed when 130 mgm./cm.? was placed each side the foil. 

The data are presented Fig. The mean value the annihilation contri- 
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Fic. absorption aluminum. 


bution the rate this geometry was 0.181 0.014 for all 
absorbers thickness over 130 This value was taken representing 
complete stoppage all positrons the absorbers and the ratios the values 
obtained for each individual position this value are expressed percentages 
the figure. lower than the value given above for point source because 
the large source area (0.5 2.0 cm.) which has marked effect the annihi- 
lation coincidence rate and very little the rate. 

The theoretical curves shown were calculated for plane sheet active 
material infinite extent and thickness covered with aluminum 
absorbers both sides. First, the path lengths from all points 
within the source the edge the absorber was computed, and the fraction 
integration over all paths was then possible compute the fraction 
homogeneous betas which would escape from the source function their 
energy. final integration over the complete positron spectrum calculated from 
the Fermi theory gave the curves shown the figure. 

seen that the experimental points are quite inconsistent with positron 
end point near 0.65 Mev. and suggest that the end point lies between 1.0 and 
1.4 Mev. 

The low specific activity the material the time the experiment and 
the restriction the geometric area the source imposed the coincidence 
arrangement made impossible use thinner source which might have given 
more clear-cut results. 

The component observed this experiment presumably due tran- 
sition from the ground state the ground state However, the 
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possibility cannot excluded that some the observed positrons are created 
nuclear the transition from the 1.856 Mev. level 
the ground state. theoretical comparison the lifetime for this 
transition (8, 27) with that the competing transitions from the 1.856 Mev. 
level (3) indicates that nuclear pair-production should not account for more than 
few per cent the positron yield. This does not conclusively eliminate this 
transition source the observed positrons since calculations this sort 
are known unreliable several orders magnitude. 

Despite the limitations this experiment, does give support the work 
McGinnis (23) who has given the value 2.07 0.2 Mev. for the 
mass difference. 

DISCUSSION 


The four experiments described above provide strong support for the decay 
scheme shown Fig. which was originally proposed the basis energy 
and intensity measurements. Although these experiments suggest that the 
1.856 Mev. level fed K-capture from the ground state the 
possibility that due transition from the spin state cannot 
ruled out. 

The log values for the transitions from the ground state the 0.0, 
0.556, 1.278, and 1.856 levels are, respectively, 6.5, 5.3, 3.4, 
and 3.4 one uses 2.07 Mev. for the mass difference. The probable 
error 0.2 the 2.07 Mev. value introduces uncertainty 0.6 into the log 
value for the 1.856 Mev. level and somewhat smaller uncertainties into the 
others. The value for the decay 4.4. This means that the transition the 
ground state roughly times less probable whiie that the two 
upper levels about times more probable than the value expected for 
allowed transition indium. This unusual situation exists also where 
the 2.25 Mev. positron decay goes 0.656 level (log 5.5) while the 
energetically favored transition the ground state missing (12). assump- 
tion that the missing transition less than 10% the intensity the 2.25 Mev. 
positron group leads log 6.5 for this transition 

has been suggested Dr. Goldhaber (11) that wide variation 
values between transitions which the basis spin and parity assignments 
are similar may expected for nuclides containing almost closed shell. The 
indium nuclei with protons all fall into this category. 

The parities all the states Fig. may expected even. This 
assignment supported the results these and other experiments and 
consistent with the proposed shell structure (12) for indium. 
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PHOTONEUTRON CROSS SECTIONS SOME 


GOLDEMBERG? LEON Katz 


ABSTRACT 


The neutron yield from photo-reactions nine elements has been measured 
function the maximum energy the University Saskatchewan 
Betatron. Cross sections are given for each the reactions calculated from 
these data. The integrated cross sections are compared with theoretical sum-rule 
predictions Levinger and Bethe and good agreement obtained with 
The data support the conclusions reached preceding paper 
Katz, and Goldemberg. The competition between and reactions 
chromium and nickel analyzed the light these experimental results. 
Good agreement with statistical theory obtained both cases. 


INTRODUCTION 


recent paper Montalbetti, Katz, and Goldemberg (21) have presented 
measurements neutron yields and cross sections for some elements, 
using method which detected the neutrons emitted these reactions. 
This work has been continued and some further measurements are here 
reported. These cover the elements lithium, chlorine, calcium, chromium, 
vanadium, nickel, bromine, lanthanum, and tantalum. Measurements 
residual activity methods some isotopes chromium are also included. 

Particular emphasis has been given nickel, since the (y, cross section 
for the element, well the cross section the isotope has 
been measured previously. Using subtraction procedure whose validity 
has been demonstrated (7), possible obtain the (y, cross section 
detailed computation the competition between neutron and 
proton emission and Ni® then gives the (y, cross section the 
natural element which compared the measured cross section. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


The apparatus used these measurements has been described 
Reference (21). The same calibration procedure was also followed the 
present measurements. 

During our early tests the sensitivity the different neutron 
energies, some the yield curves previously measured residual activity 
were checked and found have somewhat altered shape. This problem 
was not pursued further that time, but has now been more fully investigated, 
and presented appendix this paper. 

Though the full implications these changes activation curve shape 
are not understood, their effect our results (cross section curves) expected 
within the accuracy the present status the experimental techniques 
used this field. 

further indication the reproducibility activation curve shape 
obtained comparing measurements taken different laboratories the 

Manuscript received August 14, 1953. 


Contribution from the Department Physics, University Saskatchewan, Saskatoon, Sask. 
leave absence from the Departamento Fisica, Universidade Sao Paulo, Sao Paulo, 


Brazil. 


| 
_ 


50 CANADIAN JOURNAL OF PHYSICS. VOL. 32 


same element. Such comparison has already been presented for phosphorus 
(21). similar comparison the activation curve tantalum will presented 
the appendix. 

Corrections for attenuation the photon spectrum the finite sized 
samples used were applied outlined Reference (21) and the cross section 
curves were calculated using the photon difference (18) method. 


RESULTS 

The measured neutron yield curves and the cross sections are 
presented the following diagrams and the pertinent photonuclear data are 
summarized Table Results for each the elements will discussed 
separately, and for the isotopic abundance each the elements 
well their (y, thresholds will given the beginning each dis- 
cussion. thresholds were calculated using the most recent mass data. 

This the lightest element which have studied the photo-reaction 
process. The cross section presents broad peak indicating that even 
light clement, the usual picture photonuclear processes applies (8). 


NEUTRONS / 100r./ MOLE 


Li (aro) 


8 @ 6 2 6 8 20 22 24 
MAX. BETATRON ENERGY Mev. Mev. 


Neutron vield and cross section lithium. 


for the and reactions are shown Fig. Prob- 
ably large part the measured cross section due these reactions and 
reasonable separation between them and indicated the dotted lines. 
two neutrons each, the dashed curve Fig. shown for these two reactions 
needs divided two obtain their true combined cross section. 

The reaction has recently been measured residual 
this laboratory (13). Since there Coulomb barrier speak 
lithium, should expected that many protons neutrons would 
allowing for the difference thresholds. comparison the results 
shows that the and (y, reactions lithium are the same order 


magnitude. 
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Chlorine 75.4%, 13.7; 24.6%, 11.6) Fig. 2(a) and 

This measurement gives only information natural chlorine since not 
possible separate the contributions from the two isotopes. Sodium chloride 
was used the sample material, with the activation curve for chlorine being 
obtained subtraction. The reliability this procedure has been tested 
previously (7). 


MAX. BETATRON ENERGY Mev. 


Fic. Neutron and cross section chlorine. 


Calcium 96.9%, 15.8; Ca®: 0.64%, 12.7; 0.14%, 8.0; 2.1%, 
11.4) Fig. and 

The activation curve presents special features. analysis the thresholds 
the different isotopes for the (y, process together with isotopic abundance 
considerations shows that the initial bump below Mev. must attributed 
mainly and since Ca* extremely low abundance isotope. 

The cross section for the reaction has been measured 
detecting the activity the residual nuclei Summers-Gill a/. (23) and 


\ 
Fe] Summers - Gill 
wo 
2 
Wo. 


MAX. BETATRON ENERGY Mev. Mev. 


3(a) and Neutron yield and cross section calcium. The cross section obtained 
for residual activity methods indicated. 
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shown Fig. the dashed curve. The two measurements are 
excellent agreement. Subtracting this curve from our results, obtain what 
essentially the (y, cross section for the isotope. can seen that 
this isotope has peak cross section out all proportion relative 
the main isotope The abundance ratio ~1/50 whereas the peak cross 
section ratio 1/13. This discrepancy may explained the basis 
versus (y, thresholds the isotopes involved. These thresholds are given 
Table 


TABLE 


‘THRESHOLDS FOR AND REACTIONS THE 
ISOTOPES CALCIUM 


Threshold Mev.) 


Isotope 

15.8 8.1 

11.4 11.6 


Chromium 4.4%, 13.4; 83.7%, 11.54; 9.5%, 7.75; 2.4%, 
10.20) Fig. and 

The cross section for chromium presents peculiar shape owing contri- 
butions small percentage isotopes. 


nN 


NEUTRONS / 100r./ MOLE 


BETATRON 


Fic. and yield and cross section chromium. 


Fig. and are shown the activation and cross section curves for 
the reaction obtained detecting the residual activity. The 
same data the reaction also obtained residual activity 
methods are presented Fig. 6(a) and The absolute determination the 
was made comparison copper using our usual procedures (7). 

our results with statistical theory (1, and with the 
results \Waffler and Hirzel (24) shown Table All calculations refer 
16.6 Mev. excitation energy; this energy was chosen representing the 
energy the y-rays used the experiment and Hirzel. 
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(b) 


NEUTRONS MOLE 
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Fic. and (b). Activation curve and cross section the reaction 
obtained residual activity methods. 
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Fic. 6(a) and Activation curve and cross section the reaction 
obtained residual activity methods. 


TABLE 
COMPARISON EXPERIMENTS CHROMIUM WITH THEORY, 16.6 


Experiment Theory and Hirzel 


| 
| 


All thresholds for the reactions involved have been calculated using masses 
from Collins, Nier, and Johnson (4) and Duckworth (6). The theoretical 
ratio obtained and Hirzel for lower than ours, owing prob- 
ably the use wrong thresholds. 

The theoretical ratio for low factor 10. 


Vanadium 0.24%, 7.9; 99.76%, 10.8) Fig. and 
This practically 100% isotope; the activation curve and cross section 
present the usual features. 
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7(a) Neutron yield and cross section vanadium. 


Nickel 67.9%, 11.7; Ni®: 26.2%, 11.6; 1.2%, 7.8; and Ni®: 3.7% 
10.86) Fig. and 

shown Fig. the cross section for nickel very clearly indicates 
contributions from the main isotopes. Subtracting the contribution from 
measured Katz (17), the remaining part must attributed almost 
exclusively The contribution the small percentage isotopes cannot 
separated here. 


NEUTRONS / 100r./MOLE 


iv 


MAx. BETATRON ENERGY (MEV) 


Fic. and Neutron and cross section nickel. The cross section the 
reaction obtained residual activity methods indicated. 


The cross section for the emission protons has been measured for the 
natural element Halpern and Mann (9). calculated the number 
protons emitted nickel, all energies interest, assuming the validity 
the statistical theory and using the same procedure Heidmann and Bethe 
(14). The thresholds for (y, reactions used were above and for 
are those given Heidmann and Bethe (14) their paper. This computation 
was done separately for and Ni® and the results combined according 
isotopic abundance give that the natural element. 

The results are shown Fig. Halpern and Mann seems 
displaced the right about Mev.; not attribute much signifi- 
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Fic. Cross section for the reaction nickel obtained Halpern al. 
and calculated from statistical theory using the data his work. 


cance this because the results discussed the appendix. Aside from this 
displacement, the two curves are good agreement absolute value. 

The thresholds for both isotopes are quite close this case and obvious 
separation the contributions each possible. The cross section 
for the reaction has been measured residual activity methods (19) and 
plotted Fig. subtracting this from the total, get the 
cross section for the isotope 


w > 


nm 


NEUTRONS / \0Or./ MOLE 


Fic. and Neutron yield and cross section bromine. The cross section the 
(y, reaction obtained residual activity indicated. 


Lanthanum practically 106% isotope, only the characteristics 
(y, reaction quite obvious the cross section curve. Knowledge the 
cross section for the reaction would make possible separate 
out this contribution. 
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This isotope element for which the (y, reaction has been 
measured residual activity (16). Subtracting this (y, cross section from 


Neutrons U. of Tinos 
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BETATRON ENERGY Mev. 


Neutron and cross section tantalum. Data obtained residual 
and neutron detection other investigators are shown. 


the present curve the contribution from the process, starting 
obtained. These results substantiate those obtained Hanson and 


DISCUSSION RESULTS 


relevant characteristics the cross sections obtained this work are 
Table Fig. are plotted the integrated cross sections, 
together with results reported previously for number elements 
The present results reinforce the conclusions obtained previously about the 
validity the sum rules, namely that for high the emission neutrons 
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accounts for most the integrated cross section for absorption photons 
and that 0.5 seems reasonable choice for the fraction exchange- 


TABLE 
CHARACTERISTICS PHOTONUCLEAR CROSS SECTIONS 


Half-width Integrated (Mev.-barn) 


Ay, 2n) 
9.8 17.5 .004 0.020* 
5.8 19.0 .033 0.21 
5.8 17.5 105 0.60 
19.7 027 0.09 
3.6 19.0 056 0.25 
5.4 17.7 100 0.59 
6.5 16.5 050 0.34 
7.0 16.0 .175 1.3 
15.5 350 2.2 


All cross sections corrected for isotopic abundance. 
*Integrated cross section peak position. 


4 


(Mev - barn) 


INTEGRATED CROSS SECTION 


Fic. 13. Experimental values the integrated cross sections given Table Data 
contained Reference (21) are included. The solid lines represent theoretical predictions. 


force the interaction protons. For low besides the contribu- 
tion fo.dE from the process, substantial absorption photons 
must occur higher lower energies than those covered this experiment. 

The results obtained for and the comparison with the cross 
section obtained Halpern and Mann constitute direct comparison 
tween competition for neutrons and protons the same 
The agreement obtained this comparison, medium-weight nucleus, 
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for which some hope exists for the applicability the statistical model. 
quite contrary the results published and Hirzel. Further work 
the same problem with many other elements necessary before any definite 
conclusion may reached. 


APPENDIX 


During the course some experiments performed recently noticed that 
the shape the activation curve for the reaction was different 
from that obtained previously this laboratory (12, 15). All these experiments 
were made measuring the min. activity induced carbon from this reaction. 
also found that activation curves for the same elements measured 
different laboratories showed similar discrepancies shape. 

Fig. show the results for tantalum from the University 
(11) and the University Saskatchewan measurements. Measurements 
the neutron yield tantalum made the University Pennsylvania 
(10) are also included the same graph. All results have been normalized 
ours Mev. will seen that the three neutron yield curves taken 
the three laboratories are disagreement above Mev. Similarly the 
residual activity curves obtained and Saskatchewan (11, 16) are 
disagreement. However, the differences between the neutron vield curve and 
the residual activity curve plotted the lower part the diagram are 
fair agreement. similar though smaller discrepancy found the case 
when compare the measurements Illinois (5), Case Institute 
Technology (20), and Saskatchewan (16). 

would very tempting attribute the difference between the results 
taken different laboratories small changes the monitoring procedure 
the y-ray flux. This explanation is, however, partially ruled out since 
similar differences have been found between measurements taken months 
apart this laboratory using the same monitoring procedure. For this reason 
was felt that some more fundamental point was question and measure- 
ments were undertaken find what characteristics the betatron could give 
rise changes the activation curves. 

suspected that the explanation for the observed effect might result 
from changes bremsstrahlung spectrum and for this reason was felt 
that the high threshold carbon should accentuate the observed effects. 
Activation curves for the reaction were measured under different 
conditions, such changes stable orbit size and slight displacements 
the donut. 

The only noticeable effect was found when different donuts were used. 
This result not altogether unexpected, considering the rather shape 
the conventional used betatrons. shown recent article 
workers the National Bureau Standards (22) consisting 
wire mils diameter behaves far bremsstrahlung radiation 
concerned target thickness between and mils. Thus not 
unlikely that small differences between the donuts different laboratories 
and differences expander pulse have the effect producing radiation 
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different composition. should mentioned incidentally that thick target 
increases the number soft photons relative the number high-energy 
photons. This effect depresses the high-energy portion the activation curve 
since for given radiation dose there are fewer high-energy photons and thus 
less activity induced. 

should remark that the University Saskatchewan curves shown 
Fig. are the steepest all the measured curves; this would seem indicate 
that the spectrum our machine closer thin target spectrum. Moreover, 
owing the lack our knowledge the spectrum each betatron different 
times, cross sections cannot considered known better than 
absolute value; this have add the uncertainty the method 
measuring activities. 

conclusion, should emphasize that considerable care has taken 
when comparing data different laboratories, particularly differences 
between activation curves are used evaluate higher processes, such 
(y, outlined our method (7). 

Fortunately cross section shapes are not extremely sensitive changes 
activation curve shape. estimate that spite the difficulties presented 
above the peak position the cross section curve and its width 
maximum are known better than Mev. 
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ELECTRON CAPTURE SLOW SINGLY CHARGED HELIUM 
IONS 


Jackson 


ABSTRACT 


The capture electrons slow singly charged helium ions passing through 
helium gas discussed theoretically. Calculations the capture cross section 
for incident ion energies from 200 ev. 100 kev. are made using impact 
parameter method due Kohn. this energy interval the computed cross 
section found fitted very closely the empirical expression 
7.94 2.35 E(kev.). The results are compared with recent experimental 
data. The agreement between theory and experiment good, except for the 
data Keene. 


INTRODUCTION 


recent years considerable precise experimental work has been done the 
process charge exchange electron capture between atoms moving with 
relative kinetic energies from few electron volts hundreds kilovolts 
(6, 14, 17).* result, renewed interest has been aroused the theoretical 
calculation the cross sections for electron capture, with view the improve- 
ment the earlier work Brinkman and Kramers (3), Massey and Smith (12), 
and others. New calculations for electron capture protons hydrogen 
the energy range above kv. have been made independently Jackson and 
Schiff (7), and Bates and Dalgarno (1); good agreement with experiment (9, 14, 
17) was obtained. lower energies different methods calculation are neces- 
sary owing the adiabatic nature the collisions. Recently the electron cap- 
ture protons from hydrogen atoms low energies has been treated 
Dalgarno and Yadav (4), using the method ‘‘perturbed stationary 
and Kohn (10) who employed simple method. Inde- 
pendently Bates, Massey, and Stewart (2), general paper inelastic and 
charge transfer collisions low energies, have discussed impact parameter 
method similar that Kohn. They also showed that the impact parameter 
method essentially equivalent the method perturbed stationary states, 
result which was obtained independently the author (unpublished). 

The present paper concerned with the application Kohn’s impact para- 
meter method the simplest adiabatic electron capture problem for which 
there are experimental data, the capture electrons singly-charged helium 
ions passing through helium. (The calculations for the simpler problem pro- 
tons hydrogen low energies cannot compared with experiment because 
they ignore molecular effects the target material, effects which are crucial 
low energies.) Massey and Smith (12) and Dallaporta and Bonfiglioli (5) gave 
theoretical discussion the electron capture helium ions helium means 
the perturbed stationary state method. spite this earlier work was felt 
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worth while re-examine this particular problem simpler method, especi- 
ally since there were new experiments energies not covered the earlier 
calculations, and since Massey and Smith had made some rough approximations 
whose accuracy hard assess.* 


CROSS SECTION FOR RESONANCE ELECTRON CAPTURE 


For the resonance charge exchange between two like atoms (e.g. He* He) 
moving relative velocities small compared the orbital velocities the 
electrons, the cross section can written as: 


where the probability charge exchange for the two atoms passing each 
other impact parameter given equation [134] Bates, 
and Stewart (2), somewhat simplified form Kohn (10): 


where 

Here the relative velocity the two atoms, the Bohr orbit velocity 
the internuclear separation the same units, are the two lowest 
molecular potential energy curves for the appropriate molecular ion units 
(In the present case the helium molecular ion 

The differences between the present equation [2] and equation [134] Bates, 
Massey, and Stewart (2) are not significant the present problem provided the 
velocity such that (i.e. for all energies above ev.). 

characteristic result appears from equations [1] and [2]. The integral [2] 
means that oscillates very rapidly between zero and unity for small values 
computation often adequate break the range integration 
between and approximated using the average value P(p), while the 
integration beyond done numerically. Thus 


Tao 


[1] numerically. 
Two comments may made passing. The first that characteristic 


Note added proof: paper Fursov, exchange ions slow col- 
Zhur. Teoret. Fiz. 21: 1001. 1951, has recently come the author's attention. 
Fursov develops the same results Kohn, and crud. estimation the molecular potential en- 
ergy curves large nuclear separations obtains rough values for the charge exchange cross sections. 
For helium ions helium gets cross section the energy range 100 ev. kev. which 
about one and one half times that found the present paper. 
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the adiabatic resonance electron capture process that all that needed the 
molecular potential energy curves for the lowest symmetric and antisymmetric 
states the appropriate molecular ion. Once these are available the calculation 
relatively simple. The other comment that addition the resonance 
electron capture there are other, nonresonance, captures which leave the captured 
electron excited state. These nonresonance captures are small compari- 
son with the resonance capture, and will neglected subsequent comparison 
with experiment. 
CAPTURE SLOW HELIUM IONS HELIUM 


The capture electrons singly charged helium ions helium gas has been 
investigated experimentally Hasted (6) and others (11) energies below 
one kilovolt, and Keene (9) and others (11) the kilovolt region. The most 
recent experimental data are shown Fig. graph the often contradictory 
early experiments given 519 Reference (11)). The theoretical curve 
computed the basis the perturbed stationary state method Massey and 
Smith (12) and Dallaporta and Bonfiglioli (5) also shown Fig. 

Using the impact parameter method given section have calculated this 
cross section for incident ion energies from about 200 ev. 100 kev. (where 
order evaluate P(p) need know the molecular energies 
the helium molecular ion. These have been evaluated approximately 
Pauling (13) many years ago. The results given his paper are adequate for 
our purposes, especially all need are differences between the two energies. 
making use the values given Pauling, and the values 
which can obtained symmetry arguments from the atomic energy levels 
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ELECTRON CAPTURE HELIUM IONS 


the beryllium atom, was possible evaluate the integral appearing [2] 
for impact parameters from This function plotted versus 
Fig. The values for greater than about are somewhat uncertain 
owing lack knowledge and limit the cross section calculations 
energies greater than about 200 ev. For given value the values 
versus can easily found from Fig. and P(p) evaluated. 
Then the integration indicated [1] (or performed give the cross sec- 
tion the interval between 200 ev. and 100 kev. the cross section was found 
empirically fit very closely the formula, 7.94 2.35 (kev.). 

The resulting curve for the charge exchange cross section compared with 
recent experiment and earlier theoretical work Fig. The agreement with the 


CALCULATION 


0.1 10 100 
ENERGY (KEV) 


Fic. section for electron capture singly charged helium ions passing through 
helium units (0.88 cm.?) function the incident ion energy kev. 


work Hasted (6) below kev. seen good. The disagreement between 
theory and the data Keene (9) marked. could argued that the impact 
parameter method must begin fail above few thousand electron volts energy 
since the velocity approaching However, Keene’s data obviously not 
join Hasted’s data lower energies, while the theoretical curve good 
agreement there. addition, preliminary measurement 131 kev. Snitzer 
(16) gives value good enough agreement with the extra- 
polated value 3.0 from the present calculations. further argument for 
the approximate validity the present calculation the region from 
100 kev. that rough Born approximation estimate (valid above 100 kev.) 
Schiff (15) gives value 3.5 100 kev., very close the value 3.2 
obtained here. The differences between the perturbed stationary state calcula- 
tions and the present work, with its somewhat better agreement with Hasted’s 
data, can probably attributed the approximations made Massey and 
Smith (12) the course their calculations. 
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The author wishes thank Dr. Kohn for making known him the 


adiabatic impact parameter method before publication. 
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ATOMIC MASS STUDY NUCLEAR SHELL STRUCTURE 


ABSTRACT 


number new atomic masses are reported the region 126 and 
are used study nuclear stability this region. stability effects associated 
with particular neutron number are noticed but rather there broad general 
increase stability between the and 126 neutron shells. 


INTRODUCTION 


was investigated for nuclear shell effects. particular, the possible existence 
subshell effects between and nucleons was looked for, but none was 
found. Only the large effects due the wide splitting the nuclear energy 
levels coupling were clearly shown. 

The region the mass table which discussed this paper that 
which 126; and 126 neutrons represent very stable nuclear 
configurations which are attributed, the one-particle picture (10, 16, 17), 
wide splitting the nuclear energy levels and respec- 
tively. Between and 126 neutrons five subshells are filled, namely 
and These fill 92, 100, 106, 110, and 112 neutrons, 
respectively, the levels fill regularly. 

the sections follow new atomic mass values are presented, principally 
for rare earth nuclides, and nuclear stability the 126 region 
investigated. 


EXPERIMENTAL 


The new atomic mass measurements upon which this paper based were 
made with Dempster-type double focusing mass spectrograph. Both the 
instrument (3) and the method (11) have been previously described some 
detail. 

The mass differences obtained the present work are given Table 


III. COMPUTATION ATOMIC MASSES 


The data Table have been combined with existing data compile 
the list masses shown Table II, where the source each its method 
computation described footnotes. Some the rare earth elements, 
namely, europium, terbium, thulium, ytterbium, and lutetium, were not 
studied this work since they were not available metallic form. 
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TABLE 
LisT NEW ATOMIC MASS DIFFERENCES 


Mass difference 


Nuclides m.m.u. Previous measurements 

74.842 

72.344 

2As% 77.444 

2Ge” 66.9+8 

70.0+5 

2As% 

76.242 

81.0+2 

84.848 

142.9 36° 


*Some the mass differences listed this table have appeared Table Reference 12. They 
are included here since they are required for some the calculations follow. 
(4). 
Reference (8). 
Dempster, A.J. Rev. 1938. 


IV. DISCUSSION 


(a) Representation the Data 

Table contains large number isotopic masses the region 
126. The procedure for analyzing for stability effects the same that 
employed earlier paper (12). This consists plotting neutron number 
against the difference between the masses calculated from the semiempirical 
mass formula Metropolis and Reitweisner (18) and the experimental masses 
Table shown Fig. proton number plot this region gives 
new information since there are insufficient isotonic masses available. 

Since the semiempirical mass formula has compensating factors which 
allow for shell structure, should expect find large mass differences 
regions higher stability. 

second method representing the data, shown Fig. which more 
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MILLIMASS UNITS 


NEUTRON NUMBER 


Fic. plot experimental deviation from the semiempirical masses versus neutron 
number. The solid curves are the best fit for points belonging even even nuclei. The 
circles represent odd even nuclei, and the open circles represent odd odd 
nuclei. 


TABLE 
ATOMIC MASSES ADOPTED THIS PAPER 


The weighted mean for the mass the (see Reference and 
the comparison from Table The masses and ure taken from References 
and 28, respectively. 

from Table and the comparison mass Reference 12. 

Based the mass conjunction with the threshold, Reference 22, and 
the decay Reference 13. 

Derived from Table using the comparison mass Reference 

Weighted mean the value derived from Table and the value from Reference 

Weighted mean values obtained from the and the comparisons 

Table The comparison masses are taken from Reference 12. 

Weighted mean values derived from the and the comparisons 
Table using the mass Reference and the mass Reference 12. 

Derived from Table and estimated mass Rb® 84.9385 a.m.u. This estimate 

made from Fig. Reference 12. 

remeasurement old plates, Reference and the new mass Reference 12. 
Based the mass given above and the threshold, Reference 22. 
Reference 

Derived from Table and the mass The mass 94.9357 was 

obtained from the known mass Mo™, Reference 12, and the threshold, Reference 22. 

remeasurement the old plates References and 
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conventional nature, plot the binding energy per nucleon versus mass 
number. Stability effects appear this plot changes the slope the 
curve. Incidentally, the new masses reported this paper make possible 
fill the only large gap which has existed the binding energy per 
nucleon curve. 


B.E/NUCLEON MEV. 


120 140 160 180 200 
MASS NUMBER 


Fic. The binding energy per nucleon curve the region 126. 


(b) Analysis Fig. for Stability Effects 

Sharp peaks and 126 neutrons are immediately evident from Fig. 
neutron levels. addition, 126 neutrons the proton shell completed, 
and both factors contribute the stability the lead region. 

one adheres the strict one-particle model the nucleus with strong 
spin-orbit coupling, subshells are expected 92, 100, 106, 110, and 112 

sweep downward Fig. from the highly stable 
neutron configuration seems reach its lowest point the vicinity 
neutrons. slight rise appears after neutrons, but this seems part 
general rise and not associated with the filling any single subshell. 

100 neutrons the increase stability approaches its 
maximum. However, mentioned above, this gradual increase and not 
the associated with shell structure, 126 neutrons. 

106 the region 106 neutrons the stability slowly decreasing. 

110 and 112 Neutrons—At 110 and 112 neutrons the mass difference con- 
tinues decrease. This indicates that extra stability associated with 
these neutron numbers. the neutron numbers increase 126, the 
stability increases greatly, expected. 

the levels are widely split, the levels may 
fill after the group. this case neutrons would constitute subshell. 
place terms stability. 


e 
8-4 
208 
7-8 
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(c) Analysis Fig. for Stability Effects 

Fig. shows sharp discontinuities and which have previously 
been attributed (4, the filling shells these two nuclides. The effects, 
however, not show markedly Fig. since they are reduced 
division the mass number, which ranges from 130 208. This insensi- 
tive nature the binding energy per nucleon curve almost completely obscures 
the small stability effect found between the and 126 neutron shells 
Fig. 

(d) Discussion the Mid-shell Stability Effect 

The extra stability the heavier rare earth region not adequately 
explained strict one-particle picture. The smooth increase and decrease 
stability, shown Fig. not the type which characterizes the com- 
pletion nuclear shells. 

Recently, modification the strict one-particle model has been suggested 
Shalit and Goldhaber (21). this modification the structure the 
core taken into account instead regarding the core simply setting 
potential which the odd particle moves. They assume that only towards 
the closing nuclear shells exist, that is, states which one 
would expect strict one-particle model. Between shells they predict 
mixing configurations, which leads higher stability than that associated 
with the state the shell theory. 

experimental fact (25) that large electric moments exist 
between shells. These large quadrupole moments indicate highly distorted 
nucleus. The region the rare earths that region which the quadrupole 
moments are the largest all (approximately ten times the quadrupole 
moments between other shells), and hence may expect large distortion 
the nuclear core. has been pointed out Davidson and Feenberg (2) 
that distortion the core leads mixing configurations. mixed 
configurations lead higher stability, postulated Shalit and Gold- 
haber, this may account for the stability effect which have found experi- 
mentally. 

Scharff-Goldhaber (20) has plotted the first excited states 
nuclei against atomic number, and finds that between shells the first excited 
states are depressed. This depression the first excited states has also been 
qualitatively explained (21) the basis mixed configurations, which 
have pushed the first excited state lower level. Experimentally, the 
depression the first excited states shows clearly (20) the region 
126. This keeping with our present results, which indicate 
that the ground states are also depressed, although not much the first 
excited state. 

has been noted (12) that extra stability occurs between the and 
neutron shells. terms the mixing states hypothesis, should 
expect mixing much less this region since the number possible 
configurations greatly reduced. 

The lack evidence for subshells can also explained the same manner. 


its 
n- 
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This mixing configurations form state would the strict 
one-particle representation state and, hence, subshells would not exist. 
Only where the states are pure, near magic numbers, would expect shell 
structure effects. 


(e) The the Semiempirical Mass Formula 


Fig. the even isotopes the even elements are represented 
the solid curve. The odd isotopes these elements are represented 
closed circles. The position these closed circles relative the solid lines 
indicates the correctness the used Metropolis and Reitweisner 
(18) their calculations based the semiempirical mass formula. There 
random distribution the odd-A isotopes indicating that the 
used these calculations good approximation. 


(f) The Mass Difference and Possible Radioactivity 

This mass difference interest because its large magnitude, seen 
from Table 4.6 0.8 Mev. the total decay energy 
less than this figure, then should unstable. Kohman (14) 
has suggested that may radioactive but and Mulholland have 
been unable detect activity and have set half-life the event 
greater than years (19). 

possible that the decay energy the negatron branch exceeds 
Mev. (15), which case would lower than and single beta 
emission would not possible. Even this the case because the 
large mass difference, potential double-beta emitter. 


(g) Summary 

small mass effect, which represents extra stability, has been found 
between the and 126 neutron shells. This effect broad general increase 
stability which occurs the same region the recently observed increase 
the stability the first excited states. does not seem possible explain 
this phenomenon strict one-particle picture the nucleus. Reference 
made plausible explanation involving the mixing configurations 
mid-shell regions. 
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INVERTING LAPLACE TRANSFORMS THE FORM 


ABSTRACT 


Some appeal made the theory integral equations construct general 
and practical successive approximation procedure for inverting transforms 
the form given the title, where p(s) and are polynomials the degree 
p(s) being greater than the degree and h(s) any function which does 
not greatly complicate the transform. The approximations converge rapidly 
when relatively small. some cases upper and lower bounds and asymptotic 
expressions for the inverse are produced. Problems from the theory counters 
and the theory heat transfer are used illustrate the method. 


INTRODUCTION 


number physical problems can reduced the mathematical one 
inverting the Laplace transform 


less than while some simple function which need not concern 
foc the moment. 

The image U(t) such transform can usually found explicitly 
expanding powers and inverting term term. The result, however, 
rarely any practical value since usually much smaller than the values 
which one interested and there consequently large number 
terms the expansion. 

propose construct method successive approximations for the 
inversion such transform. When the method applies often leads 
natural way simple expressions for upper and lower bounds and for 
asymptotic form the desired solution. The method applies those cases 
where relatively small. Such cases are important physical applications 
because can interpreted, for example, the resolving time Geiger- 
Miiller counter the time for the passage element liquid through 
surface heat exchanger. 

Since many cases p(s) and contain only one two terms, shall 
begin illustrating the procedure with two simple applications. Then 
explain, with the aid one theorem and few remarks, how the underlying 
ideas can applied the general case. 

number authors have dealt with transforms like [1]. method proposed 
Feller (2) provides practical bounds the inverse some cases, but 
otherwise little attention seems have been paid this aspect the prob- 
lem. Hayes (4) and Sherman (8) consider such transforms when investigating 
stability problems for governed differential-difference equations. 

Manuscript received September 18, 
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Location the poles [1] and the resulting direct evaluation the inverse 
have been considered for problem econometrics James and Belz (5) 
and for problem servomechanism theory Wolfe (10). 


COUNTER THEORY 


certain problem the theory counters provides simple example 
which the method can applied. Since physical considerations motivate 
the method, shall first sketch the derivation the problem. 

assume that counter subjected random events per unit time; 
not all the events are registered because the counter’s finite resolving time 
The problem calculate the expected number registrations 
time 

Laplace transform methods are especially suited this type probability 
problem (1, and adopt Feller’s approach (2) derive expression 
for the Laplace transform 


U(t). The result (equation [6]) obtained after first consider the prob- 
ability distribution the time between succeeding registrations, and then 
the time from one registration the kth succeeding one. 

The randomness the events means: the probability that the first event, 
after any fixed instant, will occur between times and after that 
instant 


dt. 


assume the resolving time that is, assume that 
the counter registers event and only registration has taken place 
during time interval preceding the event. (When the resolving time 
the problem possesses simple solution (2).) 

the time between succeeding registrations, the distribution 
given 


t<r 
that 


the time between one registration and the kth succeeding registration 
(Fig. 1), define 


that 
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REGISTRATIONS 


Fic. Schematic diagram showing events and registrations distributed along the time 
axis. Each solid block has length and covers inactive period the counter—the one 
event shown taking place such period therefore not registered. 


from the convolution property Laplace transforms (1, 27). the 
time from any instant nof registration period, the 1)th 
succeeding registration, 


where (*) the convolution operator. Then, 


have 


Now the probability exactly registrations during time after 


Therefore 


that, from equations finally obtain 


a 


k=0 


The inverse [6] can obtained the explicit and finite form 


= == t 
Q 2 k 
co 
= 


me 
ne 
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k=1 m=0 / 
([x] the greatest integer not exceeding u(s) powers 
and inverting term term. Equation [7] fact the explicit form 
but relatively unmanageable except when small, i.e. when 
the order magnitude 

For values which one usually interested different expansion 
evidently needed. later section shall give procedure for choosing 
the expansion any particular case. For the moment sufficient note 
that goes zero with expansion powers this expression 
will therefore converge rapidly when small. The resulting series can then 
inverted term term. 

This procedure can interpreted terms the theory integral equations 
and shall see that, many cases, this interpretation leads natural 
way upper and lower bounds for the desired solution. Inverting the pro- 
posed series term term equivalent (7) solving successive sub- 
stitutions the integral equation 


0 


can easily verified. 

The first term the right side [8] simply the unperturbed solution— 
i.e. the solution for the integral term the correction for nonzero 
The successive approximations can therefore expected converge 
rapidly when the effect small. Moreover since, for and 
are positive that the integral term negative, the approximations are 
alternately upper and lower bounds for U(t); one therefore knows practice 
exactly when terminate the series. 

Expanding [6] powers and inverting term term leads 
the following new form the solution: 


k+1 [t/r] 


This series converges rapidly and its terms alternate sign. Simple but 

still useful approximations can obtained, however, directly from the 

integral equation [8] restrict our attention values for the 

approximation. straightforward matter show that this mth 

approximation 

where upper bound for even and lower one for odd 
Since usually can let and obtain immediately the 

asymptotic formula 


*The integral equation corresponding the form [7] 
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Feller (2) shcwed that this expression itself upper bound for all Exam- 


ination the bounds given equation [10] will show that, except for very 
small [11] much better approximation than 


at/(1 ar) 
which the lower bound given Feller. 


HEAT TRANSFER THEORY 


second example let consider the following typical problem the 
theory heat exchangers. Suppose that the liquid tank heated 
being passed with velocity through externally connected heat exchanger 
the surface type. Then, neglecting conductivity and transverse variations 
temperature, the temperature the liquid the exchanger 
governed 


where the (constant) temperature the heating fluid the exchanger 
and parameter depending the over-all heat transfer coefficient, the 
geometry the exchanger, and the density and specific heat the liquid. 

the liquid enters the exchanger and leaves the assump- 
tions perfect mixing the tank and losses the surroundings lead 
the boundary condition 


where the fraction tank liquid circulated per unit time. 

shall assume that the liquid initially the temperature 

For most practical purposes and are large that the first term [12] 
the Laplace transform the corresponding approximate tank temperature 


the more general case considered here obtain 


Let define 
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proceed handle u(s) the counter theory problem except that 
now add and subtract the last term the denominator and expand 
powers 

The constant taken the one real root 
[16] 


The reason for choosing this particular expansion will again follow from the 
general results the next section. For the present sufficient show 
that this expansion preserves all the desirable features our solution for 
the counter theory problem. 

Expanding indicated and inverting term term equivalent solving 
successive approximations the integral equation 


min(t. L/r) 
[17] U(t) F(t) G(u) U(t du, 
0 
where can easily found. need only: 


The upper limit [17] because, can readily verified, our 

practice L/v, the ‘passage time’ the liquid through the exchanger, and 
are both small; the integral term [17] then small correction term. 
with the counter problem, the successive approximations can therefore 
expected converge rapidly. The approximations are again alternately 
upper and lower bounds. 

The mth approximation U(t) now 


which upper bound for even and lower bound for odd. 
The required tank temperature course 


formula for the tank temperature 


Using (2) method can shown that, when the 


[14] upper bound. The bounds are interchanged when 


[19] 
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GENERAL THEORY 


each the above problems were confronted with Laplace transforms 
certain These transforms could inverted explicitly and finite 
terms, but the results were too complicated practical use. have 
used another procedure, one which generates rapidly converging successive 
approximations the solution. shall show now that the procedure can 
applied any transform the form 


where p(s) polynomial degree (we shall assume the coefficient 
1), and polynomial degree less than while h(s) some 
function which does not greatly complicate the transform. Rapid con- 
vergence the method depends being small, often the case, but 
just when small that the solution (equation [7] for the counter) 
rewrite the Laplace transform the form 


f(s) 
g(s) 


that expanding powers g(s) and inverting term term equivalent 
(7) solving the integral equation 


u(s) 


successive approximations; U(t), F(t), are course the inverses 
u(s), g(s) respectively. The essential feature our choice g(s) the 
vanishing for Properly arranged this will mean that, F(t) 
goes over into the solution the problem for the integral term 
then merely term and the successive approximations can 
expected converge rapidly for small Moreover, when, the examples 
considered earlier, know that the integral positive can conclude 
that the approximations are alternately upper and lower bounds for the 
solution. 

Our purpose now show that g(s) can always chosen that G(t) 


n 


ax) 


the denominator the given transform and then divide numerator and 
denominator this same expression can write 


p(s) 
a! 
T= 
a 
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Assuming the distinct can expand partial fractions 
obtain 


where 


the meaning that the product taken over all except 
now have 


m=1 


+ = 0, 


denominator the original transform. With appropriate changes the 
partial fractions expansion g(s), and correspondingly one can 

If, when F(t) over into the solution the problem with 
the zeros p(s) This condition uniquely determines which zeros 
must chosen. (For the counter problem can take and only 
one zero, the one the origin, chosen from the denominator 
similarly for the heat exchanger problem the only zero chosen where 
the real root [16].) Standard integral equation theory guarantees the 
convergence the approximation procedure for choice any zeros. 
However, choice other than the one required here would lead F(t) 
which could not considered satisfactory zero order approximation the 
solution, and which did not remain finite 

Finally, let consider the circumstances under which bounds for the 
solution appear the course the calculations. The integral [20] must 
positive that the integral term negative correction Whether 
not this condition satisfied given problem can decided without 
only when p(s) and are particularly simple polynomials 
When the condition not satisfied some other procedure, such that given 
Feller and referred above, must used. 

For example, interesting but not uncommon situation arises the 
counter theory problem for resolving time. Feller (2, 114) 
derives the expression 


—(s+a)r 
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for the transform the frequency function which governs the distribution 

the time between the registrations. That is, U(t)dt here the probability 

that the time between consecutive registrations between and 
obtain 


—(s+a)r —(s+a)r 
that 


The integral now negative that the integral term positive 
correction the approximation F(t). Further approximations generate only 
successively better bounds which for become 


\nt+l 


this case the asymptotic form the solution, when 


which, using Feller’s method, can shown upper bound for all 

Transforms like also arise the theory automatic controllers (6, 
and econometrics (5, 9). these applications the time constant 
usually too large for efficient use our method. 
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ABSTRACT 


The four methods calculating @p, the Debye characteristic temperature for 
specific heats, and the characteristic temperature which appears the quan- 
tum theory electrical resistance, discussed previous paper, are applied 
copper, silver, and gold. allowance has also now been made for the variation 
electrical resistance with volume take into account the effects thermal ex- 
pansion before calculated and compared with @p. Revised values for 
lithium and sodium are given. 


INTRODUCTION 


previous paper (9), which shall refer theories specific 
heat and electrical resistance metals were discussed. Four methods 
calculating the Debye characteristic temperature and the characteristic 
temperature which appears the quantum theory electrical resistance 
were compared. For convenience the given previously are reprinted 
below Tables and II, which the references equation and table num- 
bers are those paper (I). 


TABLE 
METHODS EVALUATING 


No. Method obtaining Mathematical assumption 
Cy: Equation Constant 
Equation [2] Variable parameter 
Cy/(U/T): Equation Variable parameter 
TABLE 
METHODS EVALUATING 
No. Method obtaining Mathematical assumption 
behavior 
From resistance ratio using refer- 
ence temperature, Equation Constant 
From resistance ratio using refer- 
ence temperature, Equation Variable parameter 
(dR/dT)/(R/T), Equation Constant 
(dR/dT)/(R/T), Equation Variable parameter 


VARIATION ELECTRICAL RESISTANCE WITH VOLUME 


Before and can strictly compared allowance may have 
made for the effect thermal expansion. Appendix (1) was shown 
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T 
that (as computed) very closely equal Vo) (where the 
0 


volume absolute zero), and consequently the values given apply 


essentially constant volume metal. 
Now, turning the volume expansion-coefficient metal given 


formula (7) 

AC, 

where (Griineisen’s) constant, the compressibility, and may then 


taken approximately independent temperature, and hence may 
deduced from available room temperature data. Then 


volume due solely the change with volume (cf. Mott and Jones, 
271 (15)), have 


[4] 

where the change electrical resistance due the expansion the 
metal. However, equation [4] neglects any change which may occur the 
electron-lattice interaction constant (cf. equation [5] during expan- 
sion. have sound theoretical knowledge (cf., e.g., MacDonald 
and Pearson (11)). Consequently, shall assume general 


[5] 


and determine the constant proportionality from (1) experi- 
mental measurements room temperature compressibility and change 
resistance with pressure. The calculation has also been discussed 
previously Meixner (13). 

The calculation this way has been carried through for lithium and 
sodium revise the data presented (1), and also for copper, silver, and gold. 


(i) 

contrast that most metals (cf. (1)), the electrical resistance 
lithium the volume decreases. The magnitude the change due 
the thermal expansion from the volume however, small compared 
that the other alkalies and, all, amounts only 300°K. Below 
100°K. there appreciable change the values. However, 300°K., 
the corrected values the resistance lead method 360°K. 
place was pointed out previously that the observed was 
greater than while theory predicts The correction does 


not alter this conclusion. 


(ii) Sodium 

The change the electrical resistance due thermal expansion con- 
siderably larger for sodium than for lithium and amounts some 18% 
room temperature. The values computed from the revised resistance 
data are given Fig. and differ noticeably from those given Fig. (I). 


200 


150 


100 


50 


0 100 200 300 
T eK) 


for sodium after allowance for thermal expansion. 


Methods and give which falls rapidly and becomes zero approxi- 
mately 150° and 250°K. respectively. This similar the behavior 
(see Fig. The curves using methods and are computed 
assuming that constant the region 60°K., and give value 
which increases slowly with temperature. Thus, about larger than 
the lower temperature range, but higher temperature increases 


while decreases. 
Potasstum 

Since the electrical resistance potassium shows anomalies the type 
observed rubidium above 150°K. (cf. MacDonald (10)) and the major change 
the resistance due thermal expansion takes place temperatures above 
this, have not revised the calculations for potassium. 


(iv) Copper 
indebted Dr. Dauphinee for supplying accurate measurements 
the temperature variation the electrical resistance copper. The resis- 


3 
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for copper before allowance for thermal expansion. 
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Fic. for copper after allowance for thermal expansion. 
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tance vs. temperature graph shows anomalies below the ice point. The 
resistance-rise higher temperatures however more rapid than linear with 
temperature. This behavior similar that sodium mentioned (1). 

Values computed from the measured resistance are given Fig. 
while the results the calculations after allowance has been made for the 
change resistance due thermal expansion are given Fig. The allow- 
ance for expansion the metal the case copper results which re- 
mains constant over wider temperature range. The results methods and 
are the same suggesting that the electron-lattice interaction parameter 
does remain constant the temperature interval studied. 


100 
100 200 300 400 


T (°K) 


The specific heat copper was taken from the measurements Mendels- 
sohn and Closs (14), Eucken and Werth (4), Dockerty (2), and Giauque and 
Meads (6). The variation with temperature illustrated Fig. 
the case the alkalies (see (1)) most the variation due 
the neglect terms involving the numerical evaluations. The 
agreement calculated methods and indicates that the numerical 
parameter (assumed 3R) does indeed remain constant (cf. 

The discussion given shows that the results method should 
best for direct comparison and the case copper nearly 


equal and their temperature variation similar. 


(v) 
The electrical resistance vacuum annealed strip silver was measured 


bath isopentane from room temperature down approximately 130°K 
and our Collins helium crvostat temperature. The ratio 


equals 
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Values computed from the observed resistance are given Fig. 
When allowance made for the change resistance due thermal expan- 
sion, methods and show which approximately constant the region 

Method has been applied the assumption that zero between 
150° and 200°K. The position the maximum from method shifted 
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Fic. for silver before allowance for thermal expansion. 
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Fic. for silver after allowance for thermal expansion. 
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approximately 60°K. when method applied. similar maximum ap- 
peared for lithium (cf. (1)) and was associated with change crystal 
structure. The behavior below 120°K. might possibly suggest that 
similar change occurs silver. 

Specific heat data for silver were taken from the work Keesom and Kok 
(8), Eucken, Clusius, and Woitinek (3), and Meads, Forsythe, and Giauque 


(12). The values computed using the four methods are shown Fig. 


300 
200 
=) 
100 
° 100 200 300 


comparison Figs. and shows that below slightly greater 
than However, higher temperatures contrary the predic- 
tions theory. The maximum not reflected This behavior 
similar that lithium. 


(vi) Gold 

Two specimens vacuum annealed gold wire were measured. The 
specimen with the lower residual resistance (1.8 107? the ice point 
resistance) was used for the calculation. 

Fig. shows for gold function temperature from the measured 
resistance. Since shows portion where d@/dT zero for extended 
region, methods and cannot applied. When allowance made for the 
change resistance due thermal expansion the metal the over-all varia- 

maximum occurs below 100°K. and, the case silver, might 
suggest some change structure. 

The specific heat measurements Geballe and Giauque (5) were used for 
the calculations given Fig. 10. contrast shows very little 
variation with temperature. Method was applied assuming that the results 
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for gold before allowance for thermal expansion. 
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Fic. for gold after allowance for thermal expansion. 


method give constant above 200°K., although, fact, from 
method increasing slowly throughout the range temperature. 

Since method could not applied accurate comparison can 
made between and However, appears that equal slightlv 
greater than except temperatures above 300°K. 
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100 
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Fic. for gold. 


SUMMARY 


For the highly compressible alkali metals the values calculated 
methods and are significantly altered allowance made for the varia- 
tion electrical resistance with thermal expansion. However, none the con- 
clusions given (I) changed. The values for copper, silver, and gold, 
computed after allowance has been made for thermal expansion, are more 
nearly constant than those obtained directly from the measured values the 
electrical resistance. For these three metals turns out that equal 
somewhat less than contrary the predictions theory. 
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THE EFFECT SUNRISE THE REFLECTION HEIGHT 
LOW FREQUENCY 


ABSTRACT 


has been found that when radio waves frequency kc./s. and some- 
what higher are transmitted over distances 540 km., the phase and amplitude 
the received signal change suddenly about one hour before ground sunrise 
the mid-point the transmission path. This effect results from lowering 
the nighttime reflection height km. the daytime value km. 

considering the geometry the situation, apparent that the process 
photo-ionization atmospheric molecules cannot produce the observed 
effect since the relevant wave lengths are absorbed before reaching the km. 
height. the other hand, electrons may removed from negative oxygen ions 
visible and near infrared radiation which does reach the km. level prior 
sunrise. The problem then determine the the relevant wave 
lengths the km. level, and compute the rate which electrons are 
released from negative ions. shown that the pre-sunrise effect km. 
readily explained negative ion phenomenon, and the rate fall reflection 
height discussed. 


INTRODUCTION 


has been found that when radio waves frequencies and some- 
what higher are transmitted over distances between 500 and 900 km., both 
the phase and amplitude the received signal change suddenly about one 
hour before ground sunrise when the sun’s zenith distance the mid-point 
the transmission path near 98° (1, can shown that for this position 
the sun solar ray which incident the km. level (the nighttime 
height) tangent the surface sphere about km. above the 
the positions transmitter and receiver. 


Fic. the ionosphere. 


The interpretation the phase change before sunrise that the 
height being lowered from the nighttime value km. (1). apparent 
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will shown later that photoionization atoms and molecules will 
negligible the time concerned, and other processes which can produce 
electrons must examined. has been suggested, apparently first Ratcliffe 
and later others, that the pre-sunrise effect due the release electrons 
from negative ions. what follows this problem will examined. 
necessary determine the number particles along the solar ray path the 
km. level, the residual intensity the radiation the relevant wave-length 
range after has traversed the same path, and the rate release electrons 
from negative ions. 


THE NUMBER PARTICLES ALONG THE PATH CONCERNED 


There considerable uncertainty the particle density distribution 
the earth’s upper atmosphere. Possibly the best approach adopt the 
values suggested the Panel’’ (10). The total number particles 
along the path through the atmosphere may expressed as, 


P 


where the number particles /cm.*; km.; one end the path 
the km. level and the end the path towards the sun. One may 
substitute for expression which gives the variation particle density 
with height. According Nicolet (9) 


where and are the particle densities some reference level and 
height above this level respectively, and are the scale the 
reference level and distance above the reference level respectively, and 
the scale height gradient. The path with which are concerned penetrates 
the atmosphere the km. level and proceeds upwards height km. 
Taking the former height the reference level, computing H*, Ho, and 
from the rocket data, and combining equations one may integrate the 
resulting equation give necessary divide the path into segments, 
such that each segment approximately constant. From the computed 
values appeared that five segments and hence five integrals were required. 
The same method computation has been used and outlined some detail 
Lindquist (6). check the value the path was divided into 
number small segments and the total number particles obtained 
addition. 
THE INTENSITY SOLAR RADIATION THE 
PRE-SUNRISE REFLECTION HEIGHT 

need know the number solar quanta certain wave-length 
interval incident the earth’s atmosphere, and determine the residual 
intensity this same interval after the radiation has passed along the path 
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are considering and has reached the reflection level km. From Planck’s 
law can shown that the number solar quanta the frequency interval 
entering cm.? the earth’s atmosphere per second 


where and are the solar radius and the distance from the sun the earth. 
This equation has been integrated with 6000° for quanta energy 
between and ev. and with 5000° for quanta energy between 
and ev., the temperature distribution recommended several authors 
including Bates and Massey (2). The results the integration are given 
Table expressed the number quanta/cm.?/sec. energy interval 
ev. The interval given the table chosen, since includes the wave lengths 
which are important the process the removal electrons from oxygen 
negative 


TABLE. 
SOLAR RADIATION INCIDENT THE ATMOSPHERE 


2-3 
3-4 3.2 
1.8 10% 
5-6 2.6 


The residual intensity solar radiation any depth the atmosphere 
will determined the processes absorption and scattering. apparent 
that the intensity the ultraviolet radiation which responsible for the day- 
time production the ionosphere must negligible after the radiation has 
traversed the path are considering. The total number particles along 
the path 4.6 which equivalent 1.7 km. air N.T.P. This 
thickness air reduces the intensity ultraviolet radiation near zero, 
hence the photo-ionization atoms and molecules need not considered. 

The pre-sunrise effect takes place near km. and the reflection height drops 
some km. the former level the relative importance atomic and 
molecular negative ions depends upon the particular atmospheric model 
chosen. the km. level the dominating negative ion formed from molecu- 
lar oxygen. Although the wave-length interval most effective removing 
electrons from oxygen negative ions somewhat uncertain, appears 
include the range 12000-2000A (2). However, since practically all radiation 
below 3000A absorbed along the path are considering, need not 
include radiation below this wave length. the relevant wave-length interval 
there are continuous absorption regions such are found shorter wave 
lengths owing the processes dissociation and ionization. There are discrete 
absorption bands oxygen, water vapor, and ozone, but applying Beer’s 
law and using data the total numbers the above particles and the ab- 
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sorption coefficients, can readily shown that these bands have negligible 


effect. 

must also consider the weakening the radiation due Rayleigh 
scattering. When the scattering back into the beam neglected, the intensity 
the radiation the end the path question will 


[5] 

where the scattering coefficient. have indicated above the value 
may calculated function wave length. The results are given 
Table 


TABLE 
RAYLEIGH SCATTERING FACTORS 


| 
3000 0.78 
4000 0.93 
6000 0.99 


ELECTRON DENSITIES NIGHT 


order solve this problem, data are required the rates which 
several processes proceed including attachment, detachment, and recombina- 
tion. The attachment oxygen molecules 


The value the attachment coefficient not known with certainty, but 
will adopt the figure 1.4 given Biondi (3). (Oz) 
the number densitv oxygen molecules. Following the notation Bates and 
Massey (2), the rate photodetachment/negative 


e=1 


the number incident quanta energy and the cross 
section for quanta energy convenient divide the relevant energy 
range into intervals ev. has been done Table adopt the mini- 
mum values given Bates and Massey (2), and the number inci- 
dent quanta from Tables and the rate photodetachment the order 
0.0175/cm.*/sec. There are few experimental data electron and negative 
ion densities and km. prior the pre-sunrise effect, hence must 
attempt calculate the values. The rate change electron and negative ion 
densities night may expressed the following equations, 


| 
ail 
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the nocturnal ionization rate, the values are the rate coefficients 
the various processes, and the values the densities the relevant particles. 
Since the particles involved the processes recombination and collisional 
removal electrons from negative ions are not known with certainty, have 
called these particles and The total density positive ions 
equation [8] the second and third terms are the rates removal 
electrons negative ions collision, the remaining terms the electron 
attachment and recombination rates. [9] and [10] the terms are the rates 
attachment, detachment, and recombination. Adding the three equations 


have 


view the lack concrete information the distribution and 
and consequently and the height range concerned, the second 
term the left side [11] may replaced where the total 
concentration negative ions. what follows assume that all the negative 
ions are molecular, and although this assumption will poor one the 
atmospheric model Moses and (8) correct, the error will not invalidate 
the analysis used the conclusions reached. the right side [11] 
the last two terms are replaced single term, which equivalent 
the effective recombination rate negative ions. Furthermore, replacing 
where the ratio negative ion and electron densities, and 


Integration with constant and values yields the following equation, 

The constant may determined inserting the daytime value and 
time turns out that for reasonable values and the constant 


neglected are concerned with the nighttime value when? will 
large. fact for large equation [13] becomes 


order evaluate this equation need information nighttime ionization 
rates and values expression for the may obtained from 


j 
— 
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equation [9] [10]. When state equilibrium reached, when negative 
ions are formed and destroyed the same rate, have 


Bates and Massey (2) have discussed the collisional processes which can 
electrons from negative ions. Using their data, appears that the 
nighttime values and km. are the order and 35. More 
Mitra and Jones (7) have examined the problem the nighttime 
value and have summarized the available data. The values given above 
are ¢omparable those suggested Mitra and Jones. 

There little concrete information nighttime ionization rates, and since 
the value also subject some uncertainty, equation [14] has been 
has been used, and the electronic recombination coefficient neglected. The 
results are given Table III. 


TABLE 
NIGHTTIME ELECTRON DENSITIES 


THE PRE-SUNRISE EFFECT KC./S. WAVES 


Budden (5) the electron density required reflect kc./s. 
waves over distances the range 400 4000 km. between 135 and 
the other hand somewhat higher electron density needed reflect 
waves sent vertical incidence. Since kc./s. waves transmitted over 
long distances are reflected night from height some km., can say 
the lower limit the electron density this height within the range 
values given Budden (5), and that the electron density falls less than 
values short distance below km. Referring Table noticed 
that with reasonable value for the height km., electron 
sufficient reflect kc/s. waves may maintained with relatively 
low nighttime ionization rate. 

the pre-sunrise effect due the release electrons from negative ions, 
electron density will given equation [8] with added 
term account for the photodetachment process. Then, 


this equation the detachment rate per negative ion. has been assumed 


that the dominating negative ion molecular, and electronic recombination 
which relatively slow process has been neglected. Now the sum the 


6740 1580 213 112 
9530 2240 301 158 105 
11660 2740 368 194 129 
13440 3160 425 224 149 112 
15070 3530 477 250 166 125 100 


CANADIAN JOURNAL PHYSICS. VOL. 


numbers electrons and negative ions must equal the number positive 
ions, and consider sufficiently small interval time during which the 
jonization will remain constant, then 


Substituting this the previous equation gives 


This the form 


The solution 


assigned reasonable value from Table II]. example, suppose that the 
km. and the nighttime ionization rate 15/cm.*/sec., then from Table 
and hence the order 4000/cm.* Substituting these 
figures [20] along with the value calculated previously and reasonable 
values the other quantities, can shown that the electron density would 
rise 200/cm.* with the resulting reflection kc./s. waves small time 
the order one second. This will with any reasonable range values 
the relevant parameters. 

There difficulty then explaining the lowering reflection height 
the result the release electrons from negative ions. However, must also 
consider the rate lowering. apparent from the geometry the situation 
that when the pre-sunrise effect commences, ray which incident the 
km. level has passed through considerable thickness atmosphere, and 
the value for this ray will less than the figure have used referring 
ray which incident the km. level. adding the particles along 
the path the former ray appears that the path length the order 
156 km. air N.T.P. substitute this path length and the appropriate 
Rayleigh scattering coefficients law, can seen that the intensity 
all wave lengths shorter than about 4960A effectively reduced zero. 
Under these circumstances the value much reduced and fact 
the order 0.0016. 

result the foregoing, can that when the value km. 
0.0175 and electrons are being rapidly released from negative ions, 
much smaller and electrons are being removed from negative ions more slowly 
the km. level. need know what time elapses before the value 
at-70 km. increases 0.0175, since under these circumstances have 
seen, the process electron release rapid. can shown from simple 
geometry that the zenith distance the sun must decrease 
2.5° before the ray which incident the km. level tangent sphere 
km. above the ground. Since the zenith distance the sun changes 
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per hour (1) the location where the sunrise experiments were 
performed, the time delay involved some nineteen minutes. Hence, since 
evidently the maximum period during which can vary appreciably, 
not seem that the fall reflection height during the period one hour 
due the effect discussed above. 

the foregoing analysis, was assumed that the total ionization remained 
constant during the period the sunrise effect. This may not sufficiently 
approximation. For example, the nighttime negative and positive 
some sec., then the recombination rate 
which could more than offset the ionization rate with resulting decrease 
ionization. This would have the effect reducing the rate which the 
density builds owing the release these particles from negative 

may attempt find suitable empirical equation which gives the 
required variation electron density. Such equation appears 


[2 ] Ne A ( + to)? 
where to, and are constants. Then, 


have neglected the term which expresses the rate collisional removal 
electrons from negative ions. Solving for the term and replacing 


to)? 
consider the km. level the time and take the electron and 
ion densities from Table the order 130 and 
then equation [23] vields relationship between and for some 
regarding the final value this relationship substituted 
along with the limiting values the electron density and 
then can shown that will vary appreciably over period one 


hpur. 

seems likely that the slow rate fall reflection height from 
the result least two effects, one being variable factor the 
zenith distance the sun decreases, the other reduction total ionization 
during the one-hour period the pre-sunrise phenomenon. 


REMARKS OTHER SUNRISE EFFECTS 
have discussed the sunrise effect waves propa- 
ated over long distances. When these same waves are transmitted vertical 
the amplitude begins decrease one hour before ground sunrise, 
} 


the height does not drop until some ten minutes after ground 
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sunrise (1). The explanation likely follows. The electron density produced 
ground sunrise the release electrons from negative ions not sufficient 
reflect waves transmitted vertical incidence, hence reflection 
does not take place until after sunrise. However, the number electrons pro- 
duced sufficient increase the absorption and the wave returned with 
reduced amplitude. 

Weekes and Stuart (12) have investigated sunrise phenomena waves 
near frequency 100 both vertical and oblique incidence. The 
results are generally the same those observed the lower frequency 
with apparently one difference. When the higher frequency wave 
transmitted vertical incidence, the amplitude decreases and the reflection 
height remains constant during the one hour before ground sunrise. For one 
hour after ground sunrise the height reflection falls some km.; and this 
period the amplitude remains nearly constant. possible explanation follows. 
have indicated, while electrons are being released from negative ions 
certain Jevel and the electron density building up, the amplitude 
passing wave will reduced even reflection takes place some distance 
above this level. However, once this process operates constant rate which 
may occur near ground sunrise, presumably the electron density can decrease 
again owing recombination. This effect alone would increase the amplitude 
the returned signal, but the same time the reflection height dropping 
and the amplitude decreasing, the combined effects may result nearly 
constant amplitude. 


CONCLUDING REMARKS 


From the foregoing discussion appears that the sunrise effects the 
reflection low frequency waves can explained terms the release 
electrons from negative ions. Although there are uncertainties the analysis, 
these are likely due lack concrete information nighttime ionization 
values and the rates certain reactions. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


The Effect Steady Transverse Electric Field the Viscosity Polar Liquid 
Dielectrics 


well established that the viscosity polar organic liquid which conducts appreciably 
may apparently increased the application steady electric field across the channel 
in| which the liquid flows. Andrade and Dodd (1) suggested that this increase caused 
effective reduction channel depth resulting from the formation concentrations polar 
molecules around ions the vicinity the conducting walls the channel. 

Whitehead (3) ascribed electrical effect similar cause; showed that the direct 
current between two electrodes immersed stationary transformer oil decreases over 
period about one hour steady low value. There independent evidence show that 
the electric intensity within the liquid becomes nonlinear space charges build up, 
with corresponding reduction electric current (2). 

has now been found that the fall current occurs the same rate when the liquid 
flowing between the electrodes when stationary. the current falls, the viscosity 
(which increases when the field applied) returns normal; the general nature the rate 
fall viscosity similar that the current decrease, though the two curves are 
appreciably different shape. the lowest current value, the viscosity has the same value 
that which had before the field was applied. 


CURRENT IN VISCOMETER 


RELATIVE CHANGE OF FLOW, TIME 


10 20 30 


RUN NUMBCR 


Fic. 1. Change of viscosity (full line) of amyl acetate with time for a direct potential of 600 v. across the channel 
of 0.1 mm. Successive runs are at three-minute intervals. The current is indicated by the dashed line. 


These facts are illustrated the figure, which record the change the relative time 
flow succession runs with acetate taken capillary viscometer the same 
general description that used Andrade and Dodd. 

the final asymptotic values, the removal the field for interval about three minutes 
(while viscosity measurement made) results significant change either viscosity 
current when the same field reapplied. 

These experimental results are some evidence that the space charges, indeed they exist, 

accumulate much closer the conducting walls the channel than has been suggested 
other workers. the present case, measurements the rate flow indicate that these 
ionic layers exist, they have mean depth less than 0.1 each wall, and they are not 
removed motion the liquid. also appears that the apparent increase viscosity 
strictly function the electric current, and can not spurious effect caused the 


presence space charges, 
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Corroboration tor this point view was obtained modifying Andrade and 
With the application fields different intensities, the viscosity increase 
rises approximately with the intensity the applied field, saturation value 
about and above that value there further increase this confirms 
the results Andrade and Dodd. There is, course, the time-dependent decrease these 
saturation When square wave form c.p.s. substituted for the steady field, 
however, the viscosity increase continues rise linearly with the intensity the applied 
inconsistent with the theory reduction the channel depth, which would lead 
saturation frequencies few hundred cycles per second. 
hoped carry out further experiments investigate the mechanism the viscosity 
mcrease, 
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